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ABSTRACT 



Context. Debris disks trace remnant reservoirs of leftover planetesimals in planetary systems. A handful of "warm" debris disks have 
been discovered in the last years, where emission in excess starts in the mid-infrared. An interesting subset within these warm debris 
disks are those where emission features are detected in mid-infrared spectra, which points towards the presence of //m-sized dust 
I I , grains, with temperatures above hundreds K. Given the ages of the host stars, the presence of these small grains is puzzling, and 

■ questions their origin and survival in time. 

' Aims. This study focuses on determining the mineralogy of the dust around 7 debris disks with evidence for warm dust, based on 

Spitzer/iRS spectroscopic data, in order to provide new insights into the origin of the dust grains. 

Methods. We develop and present a new radiative transfer code (Debra) dedicated to SED modeling of optically thin disks. The 
Debra code is designed such as it can determine dust composition and disk properties simultaneously. We make use of this code on 
^ ' the SEDs of seven warm debris disks, in combination with recent laboratory experiments on dust optical properties. 

5_( ' Results. We find that most, if not all, debris disks in our sample are experiencing a transient phase, suggesting a production of small 

dust grains on relatively short timescales. Dust replenishment should be efficient on timescales of months for at least three sources. 

■ From a mineralogical point of view, we find that crystalline pyroxene grains (enstatite) have small abundances compared to crystalline 
olivine grains. The main result of our study is that we find evidences for Fe-rich crystalline olivine grains (Fe / [Mg + Fe] ~ 0.2) for 
several debris disks. This finding contrasts with studies of gas-rich protoplanetary disks, where Fe-bearing crystalline grains are 
usually not observed. 

^ ■ Conclusions. The presence of Fe-rich olivine grains, and the overall differences between the mineralogy of dust in Class II disks 

' compared to debris disks suggest that the transient crystalline dust in warm debris disk is of a new generation. We discuss possible 

' crystallization routes to explain our results, and also comment on the mechanisms that may be responsible for the production of small 

CO , dust grains. 
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1 . Introduction effects, such as Poynting-Robertson drag and radiation pressure, 

that can evacuate them from the system on short timescales. This 
• --H . Several hundred stars are known to harbor debns disks The ^^^^ ^j^^ population of small dust grains is continuously 
^ . first one was discovered around Vega (|Aumannet al. || l984|) , us- replenished. Such replenishment can for example happen via de- 
^ ; ing the Infrared^onormcal Satellite (IRAS). As detailed in ^^^^^^^ collisions of km-sized bodies and will eventually result 
- . . the review by |Wyatd (l2008|), debns disks are charactenzed by ^ collisional cascade that will produce the observed excess 
photospheric emission down to mid-infrared (IR) wavelengths, g^iission, and may converge to a steady-state evolution of the 
Emission is then seen in excess at longer wavelengths, caused debris disk 
by thermal emission arising from dust, located in circumstel- 
lar belts, heated by the stellar radiation. One key characteris- Most of the known debris disks are Kuiper belt-type disk s 
tic of debris disks is that they are optically thin at all wave- (ICarpe nter etaP l2009l iLiseau et al.1 1201 Oi iLohne et al.l 1201 2h . 
lengths. A comparison with Solar System's Kuiper belt or as- showing no evidences for warm dust grains in their inner regions, 
teroid belt is often made in order to explain such systems. Since However, thanks to mid-IR observations, several unusual objects 
planetesimals and km-sized bodies do not contribute to the ther- emerged during the last years: warm debris disks. For these sys- 
mal emission, smaller dust grains that radiate more efficiendy are tems, the emission in excess starts in the mid-IR (~ 5-lOjum), 
therefore required to reproduce the observed emission in excess, suggesting the presence a material much closer to the central 
However, the smallest dust grains are also subjected to various star Our study will focus on a subset of this group: warm de- 
bris disks with emi ssion features. These very rare object s ( ~ 2% 
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G-, K-stars, whose mid-IR spectra (5-35 //m) show strong emis- 
sion features. Such emission features are usually seen in primor- 
dial disks around young Class II objects, and are associated with 
//m-sized silicate dust grains. In addition to the relatively small 
sizes of such dust grains, the presence of a strong excess at short 
wavelengths (~ 5 fim) is indicative of a high temperature for the 
dust grains, suggesting the dust is located close to the star The 
presence of these small dust grains close to the star is intriguing 
since the solution of a steady-state evolution may not be suitable 
anymore. At the ages of the systems, the smallest grains should 
have been depleted, either by radiation pressure or Poynting- 
Robertson drag. This means the disks c ould be experiencing a 
transient phase (e.g.. IWvatt et al.|[2007h . In that case, two sce- 
narios are possible, either (;) a recent collision between two or 
more bodies (planetesimals, comets, asteroids), or (//) an outer 
belt feeds the inner regions, in a very similar scenario to the Late 
Heavy Bombardment (LHB), that happened in the Solar System. 
To have a better understanding of the origin of this dust, and as- 
sess whether it is transient or not, one way is to study the Spectral 
Energy Distribution (SED) of these rare debris disks, and deter- 
mine the temperature of the dust responsible for the emission, 
and thus its distance from the star 

A powerful diagnostic to complement such study is to make 
use of the emission features that arise from (sub-) yum-sized 
grains. The width, shape and peak positions of the emission fea- 
tures prwa^e information on the composition of the dust grains 
(see iHennind l201d for a recent review). The most common 
dust species are from the silicate class, which have been de- 
tected in various environments, e.g., in the interstellar medium, 
winds of AGE stars, and objects from the Solar System. The 
main building blocks of the silicates are Si, O, Mg and Fe (then 
followed by Ca and Al). When modeling the dust mineralogy 
via spectral decomposition, one should therefore focus primar- 
ily on dust species that contained these four main ingredients. 
Silicates are an association of [Si04]''" tetrahedra, with inclu- 
sions of Mg or Fe (following cosmic abundances). Not only 
the chemical composition of the grains has an influence on the 
emission features; the internal periodic arrangement of the Si04 
tetrahedra will have a significant impact on the emission fea- 
tures. Grains with long-range order are in their crystalline form, 
as opposed to the amorphous form. As explained in detail in 
lHenniiig|(l2010l) . optical properties of crystalline and amorphous 
grains are significantly different, with multiple sharp emission 
features for the former one, while emission features are much 
broader for amorphous grains. It is therefore possible to distin- 
guish between amorphous and crystalline grains, as well as dif- 
ferent dust species. Modeling the dust mineralogy can therefore 
provide useful insights into dust composition and crystallization 
processes at stake in the debris disks, which in turn provide in- 
formation on the origin of the dust. To assess if the grains trace 
a second generation or are the remnants of the primordial, gas- 
rich disks one has to search for possible differences between the 
mineralogy of dust in Class II disks comp ared to the dust we 
observe in debris disks. For several stars, Uuhasz et al.l (|2010|) 
found pyroxene crystalline grains with a Fe fraction of about 
~ 10%, however, a result that appears to be quite general for the 
petrology of crystalline olivine grains in protop lanetary disks is 
that they are Mg-rich (and hence Fe-poor, see iBouwman et al.l 
l2008llOlofsson et al.l2010l) . The overall lack of Fe-bearing crys- 
talline grains has been studied by several authors and can be ex- 
plained by crystallization processes that take place in massive 
disks (gas-phase conde nsation, iGail 2010 or thermal annealing, 
iNuth & Johnsorill2006l) . One may wonder if we observe a sim- 
ilar result for the dust in warm debris disks. If not, we may be 



witnessing the production of a second generation of dust, where 
crystalline grains are produced by diff'erent mechanisms. 

In the following, we present in Section|2]the sample of warm 
debris disks and briefly summarize the Spitzer/iRS data reduc- 
tion. The new radiative transfer code we use is described in de- 
tails in Sect.[3] We also justify which dust species we choose to 
model the mid-IR spectra. Section|4]presents several "early" re- 
sults that we find relevant to SED modeling of debris disks, and 
we present our best fit models for individual sources in Sect.|5] 
The discussion of our results is then separated in three parts. In 
Section|6]we investigate the survival and transient nature of the 
dust grains as well as the time variability for some of the sources 
in our sample, while in Sect. [7] we emphasize on the origin of 
the crystalline grains that we detect in the dust belts. Finally, 
we briefly compare our results to what is known for objects in 
our Solar System in Sect.|8] and we summarize our findings and 
possible further improvements in Sect.|9] 

2. Spitzer/iRs data and stellar parameters 

The stellar sample (Table[Tli was gathered from the Uterature and 
investigating the Spitzer archive, looking for warm debris disks 
around solar analogs (F-, G-, K-type stars). The Spitzer/iRS spec- 
tra were retrieved f rom the Spitzer archive and reduced using the 
FEPS pipeline (see lBouwman et al 1 l2008h . 

All stellar parameters, namely T*, and t/* (eff'ective 
temperature, luminosity and distance, respectively), were taken 
from literature. To reproduce the photospheric emission we use 
Kurucz model that best matches the photometric near-IR data 
and will minimize the excess emission in the 5-8 /im range of 
the Irs data. The spectral decomposition procedure is extremely 
sensitive to the photospheric emission model used, especially for 
objects that display small mid-IR excess (e.g., HD 69830, see 
Fig.inil. 

As several Spitzer/iRS AOR (Astronomical Observation 
Request) contained both high and low spectral resolution ob- 
servations at diff'erent wavelengths (e.g. Short-Low and Long- 
High), we rebin the spectra to the minimal spectral resolution 
in order to obtain an uniform spectral resolution over the entire 
wavelength range. Doing so, the fitting process does not favor 
spectral regions with higher spectral resolution. 

Four objects in our sample are known multiple systems; 
HD 113766, HD 98800, BD+20 307 and HD 15407. For the first 
two objects, given the small separation of the different compo- 
nents and the slit size of the Irs instrument, the photospheric 
emission of the dust-free component has to be removed in order 
to consistently determine the dust location around t he star. For 
HD 1 1 3766 we use the luminosities as reported in iLisse et al.l 
(l2008l) for HD 113766 A (around which the dust is located). 
For HD 98800 we use the Keck diffract i on-lim ited observations 
of the quadruple system by iPrato et aTl (1200 ll) to model ph oto- 
spheric emission of the two components. Prato et al] (12001 1) in- 
deed provide the individual fluxes for the two spectroscopic bi- 
naries A & B (the dust being located around the spectroscopic 
binary HD 98800 B). For these two objects we will now refer to 
HD 11 3766 A and HD 98800 B. BD+20 307 is a spectroscopic 
binary (3.5 days period, Weinberger 2008) of two, almost iden- 
tical, late F-type stars and the disk is circumbinary, meaning 
no correction has to be applied to the Irs spectrum. HD 15407 
is a binary system with a 21.21" separation. The dust is lo- 
cated around the pri mary HD 1540 7 A (F5) and the secondary 
HD 15407 B (K2. ,MeHs et alJl201Cl) does not appear in the Lis 
slit, therefore no correction to the spectrum has to be performed 
for this object. 
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Table 1. Stellar sample 



Stamame 


Alias 


SpT 






d. 


Age 


AOR 


Epoch 


Ref 








[K] 


[io] 


[pc] 


[Myr] 




[d-m-y] 




HD 11 3766 A 




F4V 


6800 


4.4 


123 


10-16 


3579904 


01-03-2004 


1,2,3 


HD 69830 




KOV 


5650 


0.67 


12.6 


4000 


28830720 


14-01-2009 


4, 5,6 


BD+20 307 


HIP 8920 


GO 


5750 


1.8 


92±11 


-1000 


14416384 


15-01-2006 


7, 8,9 


HD 15407 A 




F5V 


6500 


4.5 


54.7 


80 


26122496 


09-10-2008 


10 


HD 169666 




F5 


6750 


5.16 


53.2 


2100 


15016960 


02-07-2005 


11 


HD 98800 B 


TWA4 


K5 


4600 


0.62 


44.9 


8-10 


3571968 


25-06-2004 


12, 13, 14, 15 


[GBR2007] ID 8 


2MASSJ08090250-4858172 


G 


5750 


0.8 


361 


38!^-^ 


21755136 


16-06-2007 


16, 17 



References: (1 ) IChen et al.l (l2005h. ( 2) Chen et a l] ||2006I). (3) iLisse et al.l ll2008h. (A) iBeichman etall d2005l), (5) jLisse et alj j2007l), (6) 
iBeichman etaD(l20 1l|), (7)IWeinbergelll2008.). (8).Zuckemian et al.l (l2008h. ( 9')IWeinberger et al.l(l201lh. riO')lMelis et al.l(l2010j). riDlMoor et alj 
(I2009I). (1 2) Koerne Tetal] ( l2000h . (I3) |Prato et aLllToOlD . a4) lFurlan et al.l ( 120071) . (15) ]verrier & EvansI (l2008h . (I6) lGorlova et all iimti . (17) 
iNavlor & Jelfries. (i2006h 
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Fig. 1. (2sca X A for 0.1 //m-sized grains as a function of wave- 
length for three dust species. Qsca being the scattering efficiency 
and A the albedo of the grain (see Eqn[TJ. 



3. Modeling the optically thin emission 

To model the SEDs including the Irs spectra, we use the 
"DEBris disk RAdiative transfer" code (Debra hereafter), a ra- 
diative transfer code dedicated to optically thin disks. Originally 
adapted from the "Debris d isk radiative transfer simulation tool" 
(IWolf & Hillenbrandl2005h . the code has been improved in order 
to determine the dust mineralogy as probed by the Irs spectra. 



3.1. The disk setup 

The Debra code is designed to compute the SED of a "star + 
disk" system in the optically thin regime. The star is defined 
by its Kurucz model and distance (i*, and the radiative trans- 
fer code can then compute the reprocessed (thermal and scat- 
tered) emission arising from the circumstellar dust. The disk is 
defined by the following parameters: inner and outer radii ri„ 
and rout, mass of dust Mjust and a volume density power expo- 
nent a (< 0). The density profile is assumed to be a power-law 
(Tjif) - Hif/fo)") between ri,, and rout (sampled on nr - 80 ra- 
dial points). For a given set of parameters, one can then compute 
the total emission arising fr om the disk over a large wa velength 
range (for more details, see IWolf & HillenbrandlboOSll . For all 



the sources presented in this study, we choose A^in - 0.3 jim 
and /Iniax = 100 mm. The broad wavelength range is required to 
sample the heating (stellar) and cooling (dust) radiation field to 
calculate the dust temperature correctly. 

We can question whether scattered light emission should be 
considered when trying to model mid-IR (5-35 fim) spectra. If 
the albedo A of a single particle of size s is defined as follows: 



A = 



2sca(^, S) 



2abs(^, S) + 2sca(/l, i) ' 



(1) 



where 2abs(^, ■?) and Q^cai^, s) are the absorption and scatter- 
ing efficiencies (see Sect. l3.2l for more details), then the scat- 
tered light emi ssion of this single dust gr ain is proportional to 
QsM s) X A dWolf & ffiUenbrandlllOOl . Figure[I] shows how 
this quantity decreases dramatically with increasing wavelengths 
(~ 4 orders of magnitude smaller at 5 /im than at 1 fim), for three 
diff'erent amorphous dust species, for a grain size s =0.1 yum. 
Consequently, we choose not include the contribution of scat- 
tered light in the modeling process. 

The last free parameter of our model is the grain size distri- 
bution, which we assume to be a differential power-law of in- 
dex p (dn(s) oc s''ds, p < 0) between s^in and s,^^^. One given 
model is therefore the sum of the stellar flux and thermal emis- 
sion from grains of different sizes, for several dust species, and 
is defined by four parameters ^n, rout, and p. The mass of dust 
A^dust within the disk is determined when fitting the dust com- 
position over the Irs spectral range. The fitting is performed via 
the Levenberg-Marquardt algorithm on a linear combination of 
all the thermal emission contributions from the considered dust 
species. For a single dust grain of size s, the thermal emission 
flux Femis, received by an observer at distance J*, is computed 
as follows: 



(2) 



s(r, A,s)^\ — \ X QM s)X7rx B,(A, Td), 



where Bi{A, Tj) is the Planck function at temperature T^, which 
depends on grain sizes, compositions and distance from the star r 
(see Sect. 13. 2. 2b . Finally, to determine the best fit over the entire 
SED (i.e., the disk geometry as we fl as the dust compo sition) we 
use a genetic algorithm (pikaia. ICharbonneaulll995l) over the 
four free parameters mentioned above. The code iterates over 
100 generations, with 50 individuals per generation, and fits the 
dust composition at every iteration. Each time, the goodness of 
the fit is estimated with a reduced Xr and the best fit is the one 
with the lowest 
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Table 2. Dust species used. State refers to the amorphous ("A") or crystalline ("C") form. 



Species 


State 


Chemical 
formula 


Density 
[g.cm-3] 


Mg / [Mg + Fe] 
[%] 


Shape 


Ref 


Amorphous silicate 


A 


MgFeSi04 


3.71 


50 


DHS 


1 


(olivine stoichiometry) 














Amorphous silicate 


A 


Mg2Si04 


3.88 


100 


DHS 


2 


(olivine stoichiometry) 














Amorphous silicate 


A 


MgPeSizOe 


3.20 


50 


DHS 


1 


(pyroxene stoichiometry) 














Amorphous silicate 


A 


MgSiOs 


2.71 


100 


DHS 


2 


(pyroxene stoichiometry) 














Mg-rich olivine 


C 


Mgi.85Feo.i5Si04 


3.33 


92.5 


Aerosol 


3 


(San Carlos sample) 














Fe-rich olivine 


C 


Mgi.6Feo.4Si04 


3.33 


80 


Aerosol 


3 


(Sri Lanka sample) 














Ortho-enstatite 


C 


MgSiOi 


2.78 


100 


DHS 


4 


Silica (y3-cristobalite) 


c 


SiOz 


2.27 




Aerosol 


5 


Carbon 


A 


C 


1.67 




DHS 


6 



References: (DiDorschner et al.l(ll995h . f2) |jager et aD(l2003h . (3) lTamanai & Mutschkd i l2010h . (4) lJageret aP ( 11998 j) . (5) lTamanail ( l2010h . (6) 
iJager et al.l ( ll998bl "cel600". "cel800". "cellOOO": cellulose pyrolized at 600, 800 and 1000 °C, respectively) 



3.2. The dust setup 

The Irs spectra of the objects in our sample show prominent 
emission features that are associated with warm (sub-) fim sil- 
icate dust grains. In order to reproduce the data we use sev- 
eral dust species, from amorphous (with olivine and pyrox- 
ene stoichiometries) and crystalline silicates (forsterite and en- 
statite) to crystalline y6-cristobalite silica and amorphous car- 
bon. These amorphous or crystalline dust grains, or their asso- 
ciated emission features, have been obser ved numerous times in 
disks around young Class I I objects (e.g.,lBouwman et al. 20011 



2008 



2009 



Juhasz et alJ l2010l lOhveira et al.H201 ll lOlofsson et all 



2010r ISargent et al. L2009b^, and disks around evolved 
Class III obie cts (e.g., Beichman et al. 2011. Lisse et al. 2007, 



120081 l2009i IWeinberger et all 1201 ih" "Amorphous grains of 
olivine (Mg2xFe2(i-^)Si04) and pyroxene (Mg2xFe2(i-x)Si206) 
stoichiometries are also the dominant species in the interstellar 
medium (~98%. lKemper et al.ll2005llMin et al.ll2007h . They are 
the major building blocks for the dust content in circumstellar 
environments. And yet, many uncertainti es remain on the opti- 
cal properties of these grains. For instance [Tamanai et al.l(l2006l) 
show the differences in band positions and band widths of the 
emissivity profiles for several species (crystalline forsterite and 
enstatite or amorphous Mg2Si04) when using two different lab- 
oratory techniques (KBr pellets versus aerosol measurements). 
These measured emissivity profiles are also different from ab- 
sorption efliciencies calculated using optical constant measure- 
ments. One must be aware of these unfortunate limitations, and 
interpret carefully the results from spectral decomposition when 
using several dust constituents. This is the main reason why we 
focus our study on silicate dust grains and choose not to include 
additional, more complex dust species, to limit degeneracies in 
the fitting process. Other dust components would have been in- 
cluded only if significant signatures were detected in the resid- 
uals. All the dust components used in this study are listed in 
Tabled 

3.2.1. Absorption efficiencies 

There are several ways to obtain absorption efficiencies Qabs, 
such as optical constant combined with a s cattering theory or us- 
ing laboratory extinction measurements (iHenning & Mutschkd 



I2OIOI) . When possible (and relevant, see Sect. l3.2.3l l we opt for 
optical constant as it enables us to separate the contributions 
from different grain sizes, while aerosol measurements do not 
disentangle such contributions as the sample is a mixture of 
grains with different sizes. For the amorphous silicate grains we 
use the Distribution of Hollow S pheres (DHS) scat tering theory, 
with a filling factor /^ax of 0.7. iMin et al ] (I2007h have shown 
these absorption efficiencies best match the extinction profile 
toward the galactic center at 10 and 20fim. For the crystalline 
species we also use the DHS scattering theory, but with a filling 
factor /max of 1.0, for which the absorption efficiencies usually 
best match the observed band posit ions and conipare w ell with 
other laboratory measurements (see lJuhasz et al.ll2010l) . For the 
amorphous carbon grains we use the DHS theory with a filling 
factor /max of 1. 

We compute absorption efficiencies for ten discrete grain 
sizes, logarithmically spaced between Smin = 0.1 /im and i^ax- 
For the crystalline silicate and silica grains, we limit s^ax to 
1 yum. According to crystallization models, one does not expect 
to produce pure crystals of several jim via thermal annealing. 
Grains are more likely to grow via collisional aggr egation of 
both small crystalline and amorphous material (e.g. iMin et all 
l2008h . Additionally, optical properties of large crystalline grains 
can become degenerate with optical properties of other dust 
species, since the recognizable emission features of small crys- 
tals fade away when increasing s. Limiting i,nax to 1 fj.m is in 
line with previous studies of dust in protoplanetary disks (e.g. , 
iBouwman et al. I I2OO8I lMiisz"etani20Tol loiofsson et al J 120101) 
where crystalline grains are found to be much smaller than amor- 
phous grains. Concerning the maximum grain size for amor- 
phous grains, we tested several values for i^ax and the outcomes 
of this exercise are discussed in Sect. 14. 31 



3.2.2. Temperature determination 

Young massive Class II circumstellar disks contain large amount 
of dust, but also contains gas. When modeling Irs spectra for 
these objects, one can therefore assume all the dust grains, lo- 
cated in the optically thin regions of the disks, are in thermal 
contact and have the same temperature at a given distance r. 
However debris disks are devoid of gas and the previous assump- 
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0.9 




Fig. 2. From bottom to to p: gabs values for forsterite (dashed 
grey, ISogawa et al.' 2006), aerosol measurements of forsterite 
(green, iTan ianai et al.. ,2006.). Mg -rich (blue). Fe-rich olivine 
grains (red. iTamanai & Mutschkel l2010l) and stellar subtracted 
spectrum of HD 69830 in arbitrary units (black). To ease read- 
ability, gabs values and the spectrum are scaled and offset 




8 9 10 11 12 13 

A [/ivn] 



Fig. 3. 10 /im features for photosphere subtracted spectra show- 
ing the most prominent crystalline olivine features (from bottom 
to top; HD 69830, [GBR2007] ID 8 and HD 1 13766 A). Spectra 
have been scaled and offset 



tion cannot hold anymore as no medium can efficiently support 
thermal exchange between dust grains. Consequently, we have 
to compute the temperature of the grains individually. The tem- 
perature Td of a given dust grain does not only depend on its 
distance to the star, but also on how efficiently it can absorb stel- 
lar light. The absorption efficiency being traced by Q^bsi'l, s), the 
relation between T,/ and r can be expressed as follows: 



r(s, Trf) = ^ X -I 

2 \ lKXB,i{A,Td)xQMs)dA 



(3) 



where F*(/l, T*) is the stellar surface flux at temperature T*, 
and Ri, the stellar radius. Equation|3] shows that r, and hence 
Td, depends on the absorption efficiencies used, especially in the 
near-IR where F*(/l) peaks. For the F-, G-, K-type stars of our 
sample peaks around ~ 0.5-0.6 fim, a wavelength included in 
our spectral range (/Imin = 0.3 jum). Inverting numerically Eqn.|3] 
for each dust species and for each grain sizes we are able to con- 
sistently compute their temperature depending on their distance 
to the star 



3.2.3. Forsterite and crystalline olivine grains 

Absorption efficiencies are the tools used to interpret observa- 
tions, and one should always make sure to use the most relevant 
ones. As mentioned previously, these efficiencies suffer from un- 
certainties at different levels. Laboratory measurements are un- 
fortunately subjected to various environmental effects such as 
electromagnetic polarization of the embedding mediu m (KBr, 
Csl or PE) w hen using the pellet technique (see T amanai et al.l 
120061 120091) . Additionally when using optical constant, scat- 
tering the ories have there own caveats, as discussed for ex- 
ample in iMin et all (l2007h for both Mie and DHS theories. 
By inspecting the emission features that are usually associ- 
ated with crystalline olivine grains, we can determine which 
laboratory data are best suited to reproduce the observations. 
Figure|2] shows several laboratory measurements for two wave- 



length ranges. From bottom to top, the dashed grey line shows 
absorption efficiencies calculated with the DHS theor y (with 
/max = 1) from optical constant of ISogawa et al.l (l2006l) . green, 
blue and red lines repres ent aerosol measurem ents for synthe- 
sized forsterite (Mg2Si04, Tamanai et al. 2006), Mg-rich and Fe- 
rich olivine grains (Tamanai & Mutschke 2010) . with Mg frac- 
tions (= Mg/[Mg + Fe]) of 100, 92.5 and 80 %, respectively. The 
black line is the photosphere subtracted spectrum of HD 69830. 
To ease readability, laboratory data and spectrum have been 
scaled and offset. An immediate conclusion from Figure|2] is 
that the c alculated Qabs value s of forsterite (using optical con- 
stant from lSogawa et al]l2006h do not match the observed spec- 
trum, neither around IQ fim nor at 19 and 24-jj.m. According 
to [Mutschke et al. ( 2009) some of these differences, especially 
around 10 /im, arise from a simplified treatment of the crys- 
tal anisotropy in scattering theories. At this point, it is impor- 
tant to note that HD 69830 is not the only object for which 
these differences will be a problem. Figure|3] shows the photo- 
sphere subtracted spectra (with different scalings and offsets) 
for three objects that show strong emission features associated 
with crystalline olivine grains (from bottom to top: HD 69830, 
[GBR2007] ID 8 and HD 1 13766 A). Band positions and shapes 
are overall quite similar for the three sources. For these objects 
as well, the gabs values computed with DHS will not match the 
shapes and peak positions of the observed emission features. 
Instead, the -scaled- aerosol measurements shown in Fig.|2]pro- 
vide a much better match. Consequently, we choose not to use 
the absorption efficiencies computed with the optical constant 
from ,Sogawa et al., (2006i) , and prefer to make use of the aerosol 
measurements, to represent for crystalline olivine grains in our 
modeling. The great advantage of the aerosol technique is that 
it minimizes not only the contamination by environmental ef- 
fects, but also eliminate possible computational artifacts inher- 
ited from the scattering theory used. 

The only difference between the three aerosol measurements 
shown in Figure|2]is the iron content. The Mg fractions are 100, 
92.5 and 80% for the different samples (the subscripts of Fe and 
Mg fractions of the two natural samples are determined with 
energy dispersive X-ray detection). As shown in Figure|2] the 
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Fig. 4. Peak positions of main emission features associated with crystalline olivine grains (9.8, 11, 16.2, 23-24 and 33yum), as a 
function of Fe fraction. Black circl es, red squares and green pentagons correspond to aerosol measurements of pure forsterite, Mg- 
and Fe-rich samples, respectively (iTamanai et al.ll2006t iTamanai & Mutschkell201Cih . Hatched regions correspond to the range of 
peak positions in the observed spectra of HD 69830, HD 113766 A, ID8 and HD 169666. 



iron content mostly makes a difference at wavelengths longer 
than 20 /zm, while it has very little influence around 10 /im. 
For the Fe-rich olivine grains, the 19jum feature is slightly 
broader and the feature at 23-24 jum is shifted toward longer 
wavelengths , a res ult consistent with the measurements from 
iKoike etal] ( l2003l) . To investigate this effect deeper, Figure|4] 
shows the peak positions for aerosol measurements, for several 
emission features associated with crystalline olivine grains, as a 
function of the Fe content. Peak positions for aerosol measure- 
ments for synthesized forsterite (Tamanai et al. 2006) are repre- 
sented by black ci rcles, while peak positions f or Mg- and Fe- 
rich olivine grains (ITamanai & Mutschkell20ldl) are represented 
by red squares and green pentagons, respectively. The hatched 
areas on Figure|4] show the range of peak positions for these 
specific emission features in the observed spectra of HD 69830, 
HD 1 13766 A, [GBR2007] ID 8 and HD 169666. From FigureH 
one can immediately see that the observed peak positions are 
better reproduced by Mg- or Fe-rich olivine grains, while synthe- 
sized forsterite does not seem to provide a good match (within 
a few tenths of //m). Since these three measurements are other- 
wise pretty similar, we decide to keep only the two most relevant 
samples with respect to our observations (Mg- and Fe-rich sam- 
ples), in order to limit the number of degenerate parameters in 
the modeling approach. The relevance of including crystalline 
olivine grains with a Fe fraction up to 20% is further discussed 
in Section l4~2] 

The disadvantages of using the aerosol measurements in- 
stead of optical constant are: {i) the overall level of mea- 
surements is difficult to constrain. The experiment is non- 
quantitative, in the sense that it is not possible to determine how 
many particles have entered the aerosol chamber, and thus the 
total dust mass is uncertain. And f ii), the information about the 
contribution of different grain sizes is lost. It is unfortunately not 
possible to select one finite grain size while injecting them in the 
aerosol apparatus. However, it is still possible to select a maxi- 
mum grain size, which is ~1 //m in the case of the two samples 
discussed here. One should finally be aware that extinction spec- 
tra obtained by aerosol spectroscopy contain a scattering con- 
tribution, which is however negligible since the typical particle 
sizes are much smaller compared to the wavelength of observa- 
tion. Therefore, the aerosol extinction spectra are equivalent to 
absorption in the frame of our work. 




A [/xm] 



Fig. 5. Qabs values for 0. 1 //m-sized Mg- and Fe-rich olivine 
grains (solid black line and red dashed line, respectiv ely), using 
the near-IR measurements from Zeidler et al.l (1201 ll) combin ed 
with aerosol measurements from Tamanai & Mutschkel(l2010b 



In a recent study, IZeidler et al.l (1201 ll) determined the opti- 
cal constant for several dust species in the near-IR, including 
the Mg- and Fe-rich olivines discussed in this Section. The in- 
creased Fe content (7.5 to 20%) does not dramatically change 
the absorption efficiencies in the near-IR. This is illustrated in 
Figure|5] where we show the Qabs values for Mg- and Fe-rich 
olivine 0. 1 //m-sized grains (solid black and dashed red lines, 
respectively). To compute these gabs values for different grain 
sizes, we first use the DHS theory (with /max = 1) for dif- 
ferent grain sizes (between 0.1 and 1.0 /^m) on the combina- 
tion ofbothjie^IRoptical constant from Zeidler et al. (201 1|) 
and ISogawa et al.l (l2006l) ( at longer waveleng ths). In the range 
7 < /i < 40yL(m (data from lZeidler et al.ll201 ll stop at 7 yum), we 
then replace the res ulting absorption efficiencies w ith the aerosol 
measurements from ITamanai & Mutschkel(l2010l) that are repre- 
sentative of grains smaller than 1 //m (the size distribution in the 
aerosol apparatus being unknown). Since aerosol data are non- 
quantitative, we additionally normalize their peak value to the 
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Fig. 6. Qabs values for jS-cristobalite silica grains from 
aerosol measurements (solid black line, iTamanail l2010l) . and 
0. 1 ^tm-sized grains compu ted with optical constant from 
iHenning & Mutschk3 d 19971 dashed red hne). 



peak value of the previously calculated gabs values at 11 yum, 
in order to have a more accurate overall level. The motivation 
behind this work being that Qabs values are required in the near- 
IR regime, in order to compute grain temperatures consistently 
(Eqn.El). 

3.2.4. Silica 

Silica grains have recognizable emission features around 9, 
~12.5, 16, and ~20yL(m. Several polymorphs of silica have 
been studied in detail by ISargent et al.l (12009 al) . According to 
the authors, spectroscopic observations of protoplanetary disks 
are best reproduced by annealed silica, which contains /3- 
cristobalite. They used the KBr pellet transmission measure- 
ments from iFabian et al.l (l2000l) . in which y6-cristobalite was 
formed by thermal annealing of SiOa during 5 h at a temper- 
ature of 1220K. This motivates our choice for using this par- 
ti cular polym orph to represent for silica grains. However, even 
if ISargent et al . (2009a) discussed the possible influence of the 
KBr medium before concluding it should not have altered too 
muc h the band positio ns, we prefer to use aerosol measurements 
from ITamanail (1201 Oh in order to avoid any possible contami- 
nation by the pellet medium. The maximum grain size inside 
the aerosol chamber is again 1 fim. Figure|6]shows the results of 
aerosol measurements for y6-cristobalite silica (black line). The 
major emission features peak at 9.0, 12. 5, 16.0 and 20 nm, w hich 
compares well with measurements from lFabian et al.l(l2000l 9. 1, 
12.6, 16.1 and20yum). 

To obtain gabs values, we could not follow the same pro- 
cedure as for Mg- and Fe-rich olivine grains (Sect. lXZJt since 
no data are available for wavelengths shorter than ~ 7 //m for 
y6-cristobahte grains. Theref ore, we use optical constant from 
iHenning & Mutschkd ( I1997L SiOa at 300 K), and combine them 
with optical constant from iron-free MgSiO^ in the spectral 
range 0.3-7 yum. With the DHS theory (/max = 1) we then com- 
pute 2abs values for 10 grain sizes between 0.1 and 1 fim (dashed 
red line for 0.1 j^m in Fig.|6]l. We then replace these absorption 
efliciencies with the aerosol measurements from lTamanail (|20 1 Ol) 
in the range 7-40 yum. The choice of using MgSiOs data has in 



the en d no measurable influence on the results. iKitamura et al.l 
(l2007h reviewed and studied optical properties of non-crystalline 
silica glass in different wavelength ranges, including the range 
0.1-7yum. However, in this particular spectral range, the imagi- 
nary part of the complex refractive index are extremely low {k^ ~ 
10^^). As further discussed in Sect. l4.1l such low fcj values ren- 
der the grains almost transparent to stellar light. Consequently, 
even at r < 1 AU from the star, such grains will have extremely 
low temperatures (as derived from Eqn|3) and will not contribute 
to the total emission. Therefore, silica grains need to be mixed 
with impurities, which will dominate the near-IR absorption ef- 
ficiencies. 

3.3. Determining uncertainties 

To assess the significance of our results, we compute uncertain- 
ties on relative abundances of the different dust species, as well 
as disk parameter s. This is ac hieved using the emcee package 
(Foreman-Mackev et al.ll20T2l) . which makes use of a Markov 
chain Monte Carlo (MCMC) to sample the parameter space and 
help providing uncertainties via Bayesian inference. For each of 
the 5000 runs of the MCMC procedure, values and their as- 
sociated probabilities (e^^'-^^) are stored, as well as the 4 free 
parameters for the disk geometry plus the corresponding abun- 
dances for each dust species. Probability distributions are then 
generated by projecting the sum of probabilities onto the con- 
sidered dimension (e.g., values). A gaussian fit to the distri- 
bution provides its width cr that we assume to be the uncertainty. 
An example of such a distribution is shown in Section|4|2] 



4. Early results on SED modeling 

In the course of this study, we encounter a couple of issues when 
modeling Irs spectra (e.g., poorly fitted emission features, low 
temperatures for grains with small Qabs values in the near-IR, or 
maximum grain size imax used for the amorphous dust grains). 
We detail in this Section why these issues are considered as prob- 
lematic, how we correct for them and what are the implications 
on the dust content in debris disks. In parallel, we emphasize the 
limits and possible biases of SED modeling and spectral decom- 
position for debris disks. 

4.1. "Dirty" silicates 

The fitting process used in this study highly depends on the grain 
temperatures. As explained in Sect. 13.2.21 these temperatures are 
determined using the Qabs values (from optical to mm wave- 
lengths) of the considered dust species and grain sizes. It is there- 
fore necessary to have a precise knowledge of the absorption ef- 
ficiencies. If the 2abs values are too low (e.g., ~ 10""*) in the 
optical and near-IR where the stellar emission peaks, the grains 
will become almost "transparent" to stellar light. Consequently, 
at the same distance from the star, such grains will have a much 
smaller temperature than grains with higher gabs values in the 
near-IR. If all the dust species had similar absorption efficien- 
cies in the near-IR this would simply translate into smaller dis- 
tances r to get appropriate temperatures. However, grains such 
as amorphous MgFeSi04 or MgFeSiaOe have higher gabs val- 
ues in the near-IR and are expected to be detected (as in Class II 
disks). If one consider the best fit will match crystalline features 
such as the ones from Mg- and Fe-rich olivines, in that case, the 
temperature of amorphous MgFeSi04 or MgFeSi206 grains, for 
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Fig. 7. Absorption efficiencies {s -Q.l fim) for the dust species used to model Irs spectra (except carbon) in the range 0.3-40 yum. 
Top row shows absorption efficiencies for amorphous dust grains, bottom row shows gabs values for crystalline dust grains. With 
the exception of the first two panels, dotted, solid and dashed lines represent absorption efficiencies mixed with volume fractions of 
1, 5 and 20 % carbon, respectively. 



which the Qabs values are higher in the near-IR, will most likely 
be above the sublimation temperature. 

To illustrate the problem, we model the spectrum of 
HD 69830 using the dust species listed in Tabled For Mg- and 
Fe-rich cryst a lline o livine grains, we use the near-IR data from 
IZeidler et al.l (|20I ih as described previously, which have small 
absorption efficiencies in the near-IR (~ 10 see Fig.|5]l. In 
that case, the final temperatures of 0. 1 //m-sized Mg- and Fe-rich 
crystalUne olivine grains are of about 380 and 475 K, respec- 
tively, at r = Tin - 4.4 x 10"^ AU. The corresponding tempera- 
tures for 0. 1 yum-sized amorphous MgFeSi04 grains is ~ 1970 K, 
a temperature at which they should vaporize. If one limits the rin 
value for which these grains will have a ~ 1500 K, then the 
crystalline olivine grains will be too cold and will consequently 
not contribute to the emission, resulting in a poor fit to the spec- 
trum. 

To circumvent this temperature issue, one can consider 
"dirty" silicates. In other words, one can assume that grains con- 
tain a fraction of impurities, such as transition metal ions, or are 
in contact with other dust g rains that have high ab sorption effi- 
ciencies in the near-IR (see. lOssenkopf et"anil992h . To account 
for this dirtiness, we use the Maxwell-Garnett effecti ve medium 
theory to mix different fractions of carbon ("cel600", iJager et aP 
Il998bl) with Mg- and Fe-rich crystalline olivine grains. The car- 
bon admixture was perfo rmed on the combined op tical constant 
from IZeidler et al.l (1201 II) and lSogawa et al.l(l2006l) . We then re- 
place the resulting Q^bs values with the normalized aerosol data 
for Mg- and Fe-rich olivine grains in the range 7 < /I < 40 yum 
(see last paragraph of Sect. l3.2.3] l. The choice of using carbon 
over other species (e.g., metallic iron) has very little effect on 
the results, since the primary goal is simply to heat up the grains. 
This is achieved by choosing featureless gabs values for the im- 
purities. Similar low-temperature problems will affect enstatite, 
silica (see Sect. l3.2.4l i. amorphous Mg2Si04 or MgSi03 grains 
that have low gabs values in the near-IR. For these dust species, 
from now on, we will consider that grains can be mixed with 



different volume fractions of carbon: 1, 5 and 20% to represent 
for low-, intermediate- and high-temperature grains. The need 
for dirty grains can be justified by dynamical process at stake in 
debris disks. Dust grains have most likely encountered several 
collisions during their lifetime, which renders the assumption 
of pure dust grains unlikely. Consequently we exclude from the 
modeling the absorption efficiencies with no impurities for Mg-, 
Fe-rich olivines, enstatite, silica and Fe-free amorphous grains. 
However, as we have no available diagnostics for this effect, we 
allow different levels of dirtiness for the aforementioned dust 
grains, by mixing different fractions of carbon. 

Figure|2] shows the final absorption efficiencies (for 
s =0.1 yum) used in this study to model the Irs spectra. With 
the exception of MgFeSi04 and MgFeSiiOe, the dotted, solid 
and dashed lines represent carbon admixtures of 1, 5 and 20 %, 
respectively, to account for the dirtiness of dust grains in a col- 
lisional environment. Introducing several impurities fractions 
reveals some of the unfortunate degeneracies of SED model- 
ing, especially on the disk geometry. Spectroscopic observations 
contain a finite amount of information; the "star + disk" system 
is not spatially resolved by the Irs instrument. By modeling the 
emission features, one can determine typical temperatures of the 
dust grains. However, the distance r from the star, and hence 
the geometry of the system, will depend on how efficiently these 
grains can absorb the stellar radiation. In other words, if one 
does not know the dirtiness of grains, then one should not claim 
to constrain the exact morphology of the disk. At that point, spa- 
tially resolved observations become mandatory to break down 
the number of degeneracies on the morphology of the system in 
order to refine the model. 

This underlines the overall complexity of modeling mid-IR 
spectra of debris disks in the optically thin regime. The need for 
dirty silicates, which is valid for all the sources in our sample, 
is the first result of this study. Even if the dust grains are not in 
thermal contact (one single temperature Td{r) for all species and 
sizes), pure grains as opposed to grains mixed with impurities. 
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cannot provide satisfying results, as it would result either in too 
weak contributions to the emission or, in temperatures above the 
sublimation radius for other dust species. 

4.2. Need for Fe-rich olivine grains 



From Figure|2] and the corresponding discussion (Sect. lX23T l. 
aerosol measurements of Mg- and Fe-rich olivine grains, seem 
to provide a better match to the observed Irs spectra (Fig.|3]l. 
From several previous studies about dust mineralogy in vari- 
ous circumstellar environments, we know that c rystalline olivine 
grain s are usually Mg-rich rather than Fe-rich (iBouwman et alj 
120081 lJuhasz etal.ll2010l lOlofsson et aTll2010l) . Under this as- 
sumption, we first run the Debra code including only Mg-rich 
olivine grains. This choice being motivated by the comparison of 
peak positions of several emission features (Fig.|4]i. Results for 
the two most relevant examples, HD 69830 and HD 113766 A, 
are shown in left panels of Fig. [8] For these two sources, the 
peak positions and shapes of emission features at 19 and 24jjm 
do not perfectly match the observed spectra. We therefore re- 
peat this exercise, this time including the contribution of Fe-rich 
olivine grains, and the results are shown in right panels of Fig.|8] 
Since the differences between laboratory data of Mg- and Fe- 
rich olivine grains are overall small (see Fig.|2}, the fits to the 
data do not change drastically. However, the peak positions and 
shapes of emission features at 19 and 24 jum are much better re- 
produced when including F e-rich oUvine grains (for HD 69830 
this has also been shown bv lTamanai & Mutschkell20Tol) . 

This result is further supported by the uncertainty estima- 
tion for the relative abundances. Figure|9]shows two probability 
distributions from the MCMC run for HD 113766 A. The solid 
black histogram shows the probability distribution as a function 
of the abundance of Mg-rich olivine grains, while the dashed 
red histogram shows the probability distribution as a function of 
the abundance of Fe-rich olivine grains. It immediately appears 
from Figure|9]that abundances of Fe-rich grains must be larger 
compared to abundances of Mg-rich grains in order to provide a 
good fit to the data. Following the results from this exercise, we 
decide to include aerosol measurements of Fe-rich olivine grains 
in the modeling of all the sources discussed in this study. 

As shown in both Figures|2] and |4] aerosol measurements of 
synthesized forsterite and Mg-rich olivine grains are pretty sim- 
ilar. In order to distinguish between them, one would need a pre- 
cision of a few tenths of fim (~0.1yL(m) on the peak positions 
of emission features. As explained previously, since we prefer 
to limit the number of degenerate parameters, we opt not to use 
laboratory measurements of synthesized forsterite grains (which 
seem less representative to the data, see Fig.|4]i. This choice has 
for consequence that we are not able to conclude on Fe fractions 
smaller than 7.25%. To be as conservative as possible, we will 
only claim that crystalline olivine grains have been enriched in 
Fe if the relative abundances of Fe-rich olivine grains are larger 
than those of Mg-rich olivine grains. 



4.3. Continuum emission 

When modeling the dust mineralogy of optically thick proto- 
planetary disks, a continuum component is include d in the mod- 
eling (e.g.. lJuhasz et al.l2010i lOlofsson et al.l2010l) . This contin- 
uum is originally meant to account for the emission of optically 
thick regions of the disk (the emission features arising from the 
optically thin layers). As grains larger than ~10;um do not show 
significant emission features, they are usually not included in 
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Fig. 9. ProbabiUty distributions from MCMC run for 
HD 11 3766 A, as a function of relative abundances of Mg- 
and Fe-rich olivine grains (solid black, and dashed red, respec- 
tively). 



the spectral decomposition procedure. Their actual presence in 
the proto-planetary disk is not questioned, but their contribution 
to the emission is assumed to be part of the continuum. This as- 
sumption no longer holds for debris disks that are optically thin 
at all wavelengths, therefore the influence of the maximum grain 
size has to be tested. 

First, one should note that, at a given distance r from the star, 
100 or lOOO/zm-sized grains are colder than /im-sized grains. 
The emission from large grains will therefore be more important 
at longer wavelengths (since they are colder). As grains larger 
than typically lO^um have featureless spectra, if we aim at mod- 
eling the emission features detected in spectra, the simplest as- 
sumption is to opt for 5max = 10 /im. In that case, preliminary 
tests have shown there are a couple of possible ways to com- 
pensate for the absence of large grains. The first possibility is 
an increased abundance of carbon grains, which have feature- 
less spectra as well and can serve as continuum. However, given 
that Qabs values of amorphous carbon grains are rather high in 
the near-IR, such grains will have high temperatures and will 
mostly contribute at short wavelengths. Increasing the width (or 
'"out) of the disk may circumvent this issue since the outermost 
carbon grains will be colder and can contribute as a "colder" 
continuum. However it may provide a worst fit to the emission 
features. We assume the dust composition to be homogeneous 
over the entire disk (same dust composition at all radii), there- 
fore if Tout increases by too much, some of the crystalline grains 
(e.g., Mg-, Fe-rich olivine grains) will be colder as well. This 
will boost their emission features at longer wavelengths (e.g., 
33yum). For all the objects (except HD 169666), such features 
are not detected, which practically rule out the presence of such 
cold crystalline grains. 

The second possibility to compensate for the absence of large 
grains is to mimic their presence by a shallower grain size dis- 
tribution exponent. One can understand p as a competition be- 
tween continuum emission and emission features. A steep grain 
size distribution will favor emission from small grains and their 
associated emission features, while a shallow p will favor con- 
tribution from large grains and will produce more "continuum" 
emission (large grains having featureless emission profiles). A 
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Fig. 8. Left panels: fits to stellar subtracted spectra of HD 69830 and HD 1 13766 A (top and bottom, respectively), without Fe-rich 
olivine grains. Right panels: same, as left panels including Fe-rich olivine grains in the modeling. For each panel, residuals are 
shown at the bottom (see Sect.|5]l. 



shallow p value will decrease the relative number of small pm- 
sized grains, and this will affect the total emission profile of 
amorphous dust grains. As an example, when fitting the spec- 
trum of HD 15407 A with Sniax = lOyum, we obtain a value of 
p - -1.7 0, which is far from p = -3.5 expected for a collisional 
disk (iDo hnanvi 1969). Such result would most likely require un- 
usual constraints on the dynamical processes at stake in the disk. 
We therefore estimate the most straightforward (and less biased) 
solution is to simply increase i^ax up to lOOOjum for the silicate 
and carbon amorphous grains, so that such large grains can ac- 
count for the continuum emission. Since we do not expect grain 
size distributions to be truncated for grain sizes s > 10 //m, this 
assumption appears to be a faithful description of grain size dis- 
tribution in debris disks. In FigurejlO] we show fits to stellar sub- 
tracted spectra of HD 15407 A and HD 113766 A (top and bot- 
tom panels, respectively). The solid black lines represent the flux 
contribution of silicate and carbon amorphous grains with sizes 
between 10 < i < 1000 pm. The flux fractions over the Irs range 
are of about ~60% and 38% of the total flux, for both objects. 
The conclusion of this exercise is that grains larger than 10 pm 
have to be included, when modeling disks in the optically thin 



regime. A value of iniax = 1000 pm is sufficient and provide good 
results. 

5. Individual source descriptions 

As several sources in our sample are well-studied debris disks, 
we present and discuss in the following the best fits we obtain 
for each source and compare our results to previous studies. One 
should note that we focus our study solely on Spitzer/lRS data 
and we do not include other measurements such as Mips pho- 
tometric observations. But when available, we compare the flux 
predictions of our best model with these measurements. Tables[3] 
and |4] show the disk parameters and relative abundances for the 
best models, respectively. Figure[T4] shows the fits to the stellar 
subtracted spectra as well as the residuals from the fits, which 
are computed as follows: 

DiflF [%] = 100 X ^^■°^^-^^'-^'^\ (4) 

with Fy.obs and /^v.modei being the observed and modeled fluxes, 
respectively. Figures[T5]shows the entire SEDs as modeled using 
the Debra code. We check the optical depth is always smaller 
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Table 3. Disk parameters, and corresponding uncertainties, for the best fits, as well as the maximum temperature reached 



Star 


nn 


f out 


a 


P 




T 

^ max 




[AU] 


[AU] 






[Me] 


[K] 



HD 113766 A 


0.6 ±0.1 


11.8±2.8 


-1.77±0.1 


-3.17±0.1 


7.05x10- 


-3 


1171 


HD 69830 


0.7 + 0.1 


1.5 ±0.2 


-2.47 ± 0.8 


-3.54 ±0.1 


3.70x10- 


-7 


684 


BD+20307 


0.2 ±0.1 


18.2 ±2.4 


-2.41 ±0.1 


-3.38 ±0.1 


3.75x10- 


-4 


1405 


HD 15407 A 


0.4 ±0.15 


19.2 ±5.3 


-3.00 ±0.8 


-3.12±0.1 


7.73x10- 


-5 


1531 


HD 169666 


52.7 ± 18.7 


68.3 


-2.44 ± 0.9 


-3.59 ±0.3 


1.19x10- 


-3 


215 


HD 98800 B 


1.3 ±0.3 


7.2 ± 1.9 


-1.12±0.5 


-2.99 ±01 


5.76x10- 


-2 


536 


IDS 


0.39± +0.07 


0.52 ±0.12 


-2.2 ±0.9 


-4.01 ±0.2 


2.37x10- 


-6 


864 



Table 4. Relative abundances (%) and corresponding uncertainties for the best fits 



Star 


MgFeSi04 
amorphous 


MgFeSijOfi 
amoiphous 


MgSiOj 
amorphous 


Mg2Si04 
amorphous 


Fe-rich 
olivine 


Mg-rich 
olivine 


Enstatite 


/3-cristobalite 


C 


HD 113766 A 


+ 0.6 


45.6±3.1 


+ 


0.7 


29.4 


±1.8 


7.6± 1.1 


1.7 = 


tO.9 


0.8 + 0.8 


+ 


0.7 


14.9 ±2.6 


HD 69830 


24.3 ±5.9 


29.9 ±0.7 


+ 


0.7 


3.7 


+4.2 
-3.7 


6.5 ±1.8 


1.4: 


tO.6 


2.1 ±0.7 


2.6 d 


:2.6 


29.5 ±6.2 


BD+20307 


68.6 ±2.5 


15.4 ±0.9 


+ 


0.6 


7.3 I 


1 1.7 


4.0 ±0.7 


0.6 


+0.7 
-0.6 


4.1 ±0.8 


+ 


0.7 


+ 0.7 


HD 15407 A 


+ 0.6 


25.0 ±5.4 


28.0 


±7.2 


+ 


0.6 


+ 0.5 


+ 


0.8 


0±0.6 


17.3 


±2.1 


29.7 ± 12.8 


HD 169666 


+ 0.6 


+ 0.7 


+ 


0.7 


65.1 


±0.7 


+ 0.8 


6.9 = 


t3.4 


7.7 ± 1.2 


3.6 d 


:0.7 


16.7 ±6.5 


HD 98800 B 


25.9 ±0.6 


28.0 ±0.7 


+ 


0.5 


22.6 


±6.7 


+ 0.8 


0.2 


+0.8 


5.2 ±1.5 


+ 


0.7 


18.2 ±6.5 


ID 8 


26.4 ± 12.0 


+ 0.6 


+ 


0.7 


60.8 ± 13.2 


2.0± 1.1 


0.3 


-M 

-0 3 


+ 0.5 


+ 


0.6 


10.5 ±3.1 




Fig. 10. Fits to stellar subtracted spectra of HD 15407 A (top) 
and HD 1 13766 A (bottom). Fits are in red, data in dashed black. 
Contribution of grains larger than 10 fim is shown in solid black. 



than unity meaning all sources can be modeled within the optical 
thin regime. 

5.1. HD 11 3766 A 

Previous studies: The Irs spectrum of HD 11 3766 A (the pri- 
mary of the binary system) shows several strong emission fea- 
tures attributed to crystalline grains (9.95, 11.0, 16.25 and 
23.5 yum). The circump rimary dust orbitin g the F3-type star has 
been studied in detail bv lLisse et al.l (l2008h . To model the Spitzer 
data, the authors used several species among a vast library of 80 
different chemical components. Among them, they found two 
types of lV[g2Si04 forsterite, labeled as "ForsteriteKoike" and 
"Forsterite038". Unfortunately, the authors do not discuss how 



different these measurements are (in Lisse et al.ll2008l or in the 
supporting online material of iLisse et al. 1 12006b . or the implica- 
tions of why both species are required. In addition to olivine and 
pyroxene dust species, the authors state they need additional dust 
species to model the spectrum, such as ferrosilite (Fe2Si206), 
diopside (CaMgSiaOe), phyllosilicates (smectite notronite), wa- 
ter ice and metal sulfides (ningerite MgioFegoS). Concerning the 
latter dust species, to our knowledge, there is a major gap in the 
wavelength range of laboratory measurements. Opt ical constan t 
are available between 0.3 and 1.1 um (iEgan & Hilgemanlll977l) 
and between 10 and 500 yum jBegemann et al. 19941) . The Q^bs 
values remain unknown both in the near-IR (rendering the tem- 
perature determination uncertain), and in the 5-10yL(m range, i.e., 
in the Irs spectral range. Additionally, in the range 10-40yL(m 
only one emission feature is visible and is very broad (~26- 
35 fim). We estimate that the lack of a clear, narrow emission 
feature makes the identification of this particular dust species 
difficult, the risk being to use metal sulfides as a proxy for con- 
tinuum emission (se e Sect. 14.311. Finall y, concerning the geom- 
etry of the dust belt, iLisse et al] ( l2008h found the dust to be lo- 
cated at 1 .8 AU from the star, with a width of 0.4 AU with typical 
temperatures between 420 and 490 K. The grain size distribution 
for their best fit has an exponent of -3.5+0.2 for grain sizes be- 
tween 0.1 and 20 /im, and finally the authors found a mass of 
~ 5 X 10"^ Me for grain sizes up to 20 //m. 

The amorphous phase: According to our model, the abun- 
dances of amorphous grains with olivine and pyroxene stoi- 
chiometrics are of the same order of magnitude (46 versus 29%, 
respectively), with an edge for gr ains of the pyroxen e group, a re- 
sult different from the analysis of lLisse et al ] (120081) who found a 
strong depletion of amorphous grains with pyroxene stoichiom- 
etry. We find that amorphous grains of the pyroxene group are 
all Fe-bearing while grains from the olivine group are Fe-free 
(no ]VIgFeSi04). The abundance of amorphous carbon grains is 
of about 15%. Grains between 10 and 1000/im in size contribute 
to about 38% of the total emission over the entire Irs range. 

The crystalline phase: According to our model, the crys- 
talline content simply consists of Fe-rich olivine grains (~8%) 
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and a small amount of Mg-rich olivine grains (~2%), both 
species being mixed with 20% of carbon. Enstatite is tentatively 
detected at 0.8% (±0.8). No silica grains are used (and needed) 
in the best fit model. Inspecting the residuals of the best fit, we 
do not see any con clusive evid ences for additional species such 
as the ones used bv lLisse et al.l (|2008i). 

Disk geometry: The dust belt is found to be rather extended, 
with Tin =0.6 AU and rout = 11.8 AU. The grain size distribution 
follows a power-law with an exponent of /:> = -3.17 and we 
obtain a total mass of dust of ~ 7 x lO"-' M© for grain sizes up to 
lOOO^/m. 

Additional obs ervations and comments: In a recent study, 
I Smith et aP (1201 2h presented VLTI/Midi 10 //m observations of 
HD 11 3766 A. Modeling the observed visibilities, the authors 
find the debris disk to be partially resolved and their observa- 
tions are consistent with a Gaussian profile with a FWHM of 
about 1.2-1. 6 AU (hence, r~0.7AU). Unfortunately, since the 
interferometric observations are partially resolved, the exact lo- 
cation of the inner radius r\n can not be further constrained. Visir 
images from lSmith et al.l(l2012l) additionally place an upper limit 
on the extension of the disk at 17 AU from the central star, a re- 
sult consistent with our best model (rout = 12 AU). 

Our best fit returns a 70yum flux of 269 mJy, to b e com- 
pared to the Mips flux of 390+28 mJy dChen et alj|201 U re vised 
from 350+35 mJy originally reported in IChen et al. I l2005h . The 
~ 120mJy difference can be explained in two different ways. 
First, one may consider time variability between Irs and Mips 
observations. In the case of HD 11 3766 A, this explanation is 
quite unlikely as only a few days separate the two observations 
(23"' of February 2004 and I ''' of March 2004 for Mips and Irs, 
respectively). Seco nd, an outer dust b elt might surround the in- 
ner warm dust belt. iLisse et al ] (120081) considered an equivalent 
scenario as their model suggests the presence of two additional 
dust belts, one as far out as 9 AU where water ice may be located 
and another one at 30-80 AU. We note that according to our best 
model, the dust belt is rather extended (up to 12 AU), suggesting 
the need for cold dust grains, which may support the possibility 
of an outer belt. 



5.2. HD 69830 

Previous studies: In addition of showing emission features 
in its Irs spectrum, HD 69830 is known to host 3 Neptune- 
sized planets (Lo vis et al.ii2006i) . with semi-major axis of -0.08, 
0.186 and 0.63 AU, with respective minimum masses of 10.2, 
11.8 and 18.1 M©. The mineralogy of th e dust content has 
been studied in detail by iLisse et aP (l2007h . who modeled the 
Irs spectrum in a similar way as the one of HD 11 3766 A. 
According t o their model, the authors detected several species. 
Similarly to lLisse et"an (l2008h . they used two types of Mg2Si04 
forsterite ("ForsteriteKoike" and "Forsterite038"). Fayalite, the 
Fe-rich end member of the olivine group was also included 
in their model. No amorphous grains with a pyroxene stoi- 
chiometry were detected, but bronzite (Mgi_iFej:Si206) fer- 
rosilite (FeaSiaOg) and diopside (CaMgSiaOe) were necessary 
for their model to reproduce the data. Finally, the authors state 
they detected several carbonate species, including magnesite 
(Mg COa), dolom i te (Ca MgCiOe), water ice and amorphous car- 
bon. iLisse et al ] (I2007h concluded the dust belt is located at 
1 .04+0.1 2 AU from the central s tar. According to the stability 
analysis from iLovis et aP (l2006l) . this value is consistent with 
the 2:1 and 5:2 orbital resonances with the outer Neptune-sized 
planet. The typical temperature for the dust is ~ 400 K and the 
grain size distribution exponent for grain sizes between 0. 1 and 



10/im is -3.9+0.2, with a dust mass for these grain sizes of 
~ 5 X 10"** Me. 

In a recent study, iBeichman et al.l (1201 ih investigated the 
time variability of the dust around HD 69830. With new spec- 
trosco pic observations and using the model from iLisse et all 
(l2007h . they led a similar mineralogical analysis. Their model 
confirmed the previous detections of amorphous grains with 
olivine stoichiometry, both species of "ForsteriteKoike" and 
"Forsterite038" and fayalite. For the pyroxene group, only 
bronzite was detected in the new analysis (while ferrosilite and 
diopside went from "detected" to "marginally detected" in the 
new analysis). The carbonates species of magnesite and dolomite 
were also conside red as marg i nally detected, as opposed to the 
best fit results of iLisse et alT(l2007l) . Sulfides (FegoMgioS) for 
which an upper limit were found in the 2007 analysis was de- 
tected in the latest best fit model. According to Beichma n et al.l 
(|20TT1) . the most significant difference between the two mod- 
els, is that water ice was no longer needed to reproduce the 
spectrum, therefore concluding the debris disk to be really dry. 
Based on the main mineralogical markers from both analysis, 
the authors state that the dust has not evolved significantly in 
time, but instead that the aforementioned differences arise from 
higher signal-to-noise ratio (SNR hereafter) of the new obser- 
vations and a better data calibration for the Irs spectra. This 
comforts us in our conclusion that spectral decomposition is a 
non-trivial problem, strongly dependent on several aspects, from 
SNR in the data, to imperfect knowledge of absorption efficien- 
cies. This is why we prefer to use less dust species but focus 
on those that have been intensively studied in laboratory exper- 
im ents, and have a go od coverage of their optical data. Finally, 
iBeichman et a mentioned that no gas lines could be con- 

vincingly detected in their data. 

The amorphous phase: According to our model, the amor- 
phous dust content is best described by Fe-bearing amorphous 
grains with olivine stoichiometry (24%) and Fe-bearing amor- 
phous grains with pyroxene stoichiometry (30%). No MgSi03 
amorphous grains are required in the best fit model, while a small 
fraction (~ 4%) of Mg2Si04 grains are found. The abundance of 
amorphous carbon is of about 30%. Grains with sizes between 
10 and 1000/^m amounts for 15% of the total flux in the Irs 
range. 

The crystalline phase: The total crystallinity fraction of 
~13% is the sum of contributions from Fe-rich olivine grains 
(~7%, mixed with 20% carbon), Mg-rich olivine grains (~1%, 
mixed with 1% carbon), enstatite (2%, mixed with 5% car- 
bon), and yS-cristobalite (3%, most of it being mixed with 1% 
carbon). Inspecting the residuals between the modeled and ob- 
served fluxes, we do not see evidences for other emission fea- 
tures. 

Disk geometry: The dust belt is located between 0.7 and 
1.5 AU, values that are not in contradiction with the position of 
the outermost Neptune-sized planet. Since we assume the disk 
to be circular (eccentricity e - 0), the so lution we fi n d is in 
agreement with the stability analysis from iLovis et aH (l2006l) . 
The grain size distribution exponent is p = -3.54, with a total 
dust mass Mjust = 3.7 x 10"^ M@ for sizes up to 1000 /.m. 

Ad ditional observations and comments: Beichma n et al.l 
(l2005h reported a Spitzer/Mips flux of 19+3 mJy at lOjjm. Our 
model returns a slightly higher flux of 22 mJy, within the un- 
certainties of the lOjjm observations. This suggests there is no 
colder material that contribute up to 70 /im and that the system 
can be successfully described with one single dust belt close 
to the star. Our best fit model is consi stent with the VLT I/Midi 
observation and VLT/Visir imaging of ISmith et al.l (l2009l) . who 
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constrain the warm dust in HD 69830 system to lie within 0.05- 
2.4 AU. 

5.3. BD+20307 

Previous studies: The dust in the debris disk of BD+20 307 sur- 
rounds a s pectroscopic binary of similar F-type s tars, a s ex- 
plained by IWeinbergeJ (l2008l) . IZuckerman et al.l (l2008l) esti- 
mated an age of a few Gyr for this binary system. The Irs 
spectrum shows prominent feature at lOyum, with a shoulder 
at 11 yum. At longer wavelengths, smaller emission features are 
visible at 16.5, 18.7 a nd 23.5-24 ^um. No dete ctable feature 
is seen around 33 jum. IWeinberger et al.l (1201 ll) modeled the 
spectrum of this object, using amorphous grains with olivine 
(MgFeSi04) and pyroxene (MgFeSiaOg) stoichiometrics plus 
forsterite grains. According to their modeling, the dust is lo- 
cated at 0.85 AU from the binary, with typical temperatures be- 
tween 450 and 500 K for the dust grains. The inferred crys- 
tallinity fraction is ~ 20% in forsterite grains. Finally, their 
model slightly over-predicts the Mips flux at 70 /im (measured: 
27.8 mJy, against 40.2 mJy for their model). 

The amorphous phase: The amorphous dust content is dom- 
inated by Fe-bearing grains from the olivine group (69%), but 
also consists of Fe-free amorphous grains with olivine stoi- 
chiometry (7%) and Fe-bearing grains with pyroxene stoichiom- 
etry (15%). No amorphous carbon grains are required in the best 
fit model. Grains with sizes between 10 and 1000/im contribute 
to 20% of the total flux in the 5-35 /im range. 

The crystalline phase: The remaining relative abundances 
are divided between Fe-rich and Mg-rich olivine grains (4 and 
0.6%, respectively) both with a 20% mixed carbon fraction and 
warm enstatite (4%) in contact with 20% of carbon. 

Disk geometry: We find the dust belt to be located between 
Tin =0.2 AU and rom =18 AU. According to Weinberger (2008), 
the spectroscopic binary has a separation of 0.06 AU at present 
days, meaning our estimation of ri„ is compatible with the ge- 
ometry of the system. The grain size distribution exponent p is 
-3.38 and the total dust mass is of about 3.7 x lO^^M® for sizes 
up to lOOOyLfm. 

Additional observations and comments: Our model slightly 
over-predicts the Mips observation, as we obtain a lOfim flux 
of 36mJy, to be compared with 27.8 mJy. Interestingly, our fit 
to the lO/um fea t ure do es not match perfectly the data, as in 
IWeinberger et al.l (l201 ih . The feature at lO/jm looks "primitive" 
(or pristine, quite similar to ISM dust), with a shoulder around 
1 1 fim but no clear emission features on top of it, while our best 
fit model shows small features caused by the presence of crys- 
talline olivine grains and enstatite. The inferred dust composition 
contains Fe-rich olivine grains because of the features detected at 
longer wavelengths (16.5 and 24 //m). In other words, the spec- 
trum seems to show two different dust compositions, a first one 
with very low crystallinity and a second one with an increased 
crystallinity fraction. To reproduce the strong, almost pristine 
lOyum feature, amorphous grains need to be warm enough to re- 
produce the data, therefore they need to be lo cated close to the 
star (0 .2 AU in our case, 0.85 AU according to IWeinberger et all 
As we assume the disk to be homogenous, as soon as 
crystalline olivine grains are present they will contribute around 
10 //m. One may argue that the temperature of the olivine grains 
could in principle be lower and thus contribute only at longer 
wavelengths. In our best fit, the mixed carbon fraction for the 
Fe-rich olivine grains is 20%, while we allow it to be 1 or 5% 
(which would provide smaller temperatures). Such solution is 
apparently not providing a better fit to the data. Allowing smaller 



fractions of carbon would basically mean that olivine grains are 
in their pure form, an unlikely prediction for dust grains that 
have most likely encountered collisions. The small uncertainties 
on the relative abundances of Fe-rich grains confirm that the fit 
cannot be improved by much, with the assumption of one sin- 
gle belt. This overall means we should consider alternative solu- 
tions. 

As discussed previously, spectroscopic observations are un- 
resolved and do not contain spatial informations. The simplest 
solution is always to assume that the emission probed by Irs is 
arising from one single dust belt. With this assumption, it ap- 
pears the fit to the data could be improved. Additionally, the 
shape of the emission features, and thus the mineralogy suggests 
a two-components dust composition. Even though the 70 //m 
Mips observations do not suggest the presence of colder mate- 
rial, one can imagine the emission in excess is arising from two 
separated dust belts: a warm inner dust belt, containing mainly 
amorphous grains and a second, slightly colder dust belt with a 
higher crystallinity fraction. This possible scenario is also sup- 
ported by the extension of the dust belt, up to -18 AU, which 
suggests the presence of rather cold dust, further away from the 
star. 



5.4. HD 15407 A 

HD 15407 A is the primary of a binary sys tem (21.25" sepa - 
ration), with an age estimate of ~80Myr (iMelis et alJl2010l) . 
The Irs spectrum mostly displays one strong emission feature 
peaking at ~ 9.1 pm that can be attributed to y6-cristobalite Si02 
grains. 

The amorphous phase: No amorphous grains with olivine 
stoichiometry are required in order to reproduce the spectrum of 
HD 15407 A. Fe-bearing amorphous grains from the pyroxene 
group amount for 25%, and Fe-free grains represent 28% of the 
total abundance. The abundance of amorphous carbon is 29.7%. 
Grains with sizes larger than 10 pm contribute significantly to 
the total emission (60%, see Fig.fTOli. 

The crystalline phase: Beside the aforementioned species, 
the only remaining dust composition required to reproduce the 
data is /3-cristobalite, for which we find an abundance of 17% 
with a mixed carbon fraction of 1%. Even though the emission 
feature at 12.5 pm is slightly under-estimated (but still visible in 
the model), the 16 pm feature is well -reproduced, as well as the 
broa der feature at 20 jUm. As shown in Fig.|6](as well as Fig. 2 to 
4 of lSargent et al.ll2b09ab . yS-cristobalite (annealed silica) is one 
of the few polymorphs of Si02 that display an emission feature 
at 16 pm (coesite also displays a similar emission feature at these 
wavelengths). Therefore, it is difficult to firmly conclude on the 
exact nature of the Si02 polymorph we detect in the disk around 
HD 15407 A, but we can narrow down possible polymorphs ei- 
ther to j6-cristobalite or coesite. No additional crystalline grains 
are found in the best fit, and the residuals do not show evidences 
for other emission features. 

Disk geometry: For the best model, the dust is located in an 
extended belt, between 0.4 and 19.2 AU. The grain size distribu- 
tion follows a power-law distribution with a slope in p - -3.1, 
and the total dust mass is ~ 7.7 x 10"^ Mg, for grain with sizes 
smaller than lOOO/^m. 

Additional observations and comments: The large extension 
of the dust belt (up to ~19 AU) underlines the need for cold dust. 
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Fig. 11. Fits to stellar subtracted spectra of HD 169666 (fits in red, data in dashed black), for the "cold" (left panel, Tmax ~ 215 K, 
r ~ 53 AU) and "warm" solutions (right panel, Tmax ~ 380 K, r ~ 13 AU). 



5.5. HD 169666 

The emission in excess seen in the SED of HD 169666 starts 
at longer wavelengths (~ lOyum) compared to the other objects 
of our sample, but still present emission features characteristic 
of small grains, especially at 24, 28 and 33 yum , that c an be at- 
tributed to crystalline olivine grains. iMoor et al.l(l2009l) reported 
a Mips flux of 21.8+4.0mJy at 70 //m. 

The amorphous phase: The composition of the amorphous 
content is solely represented by Fe-free grains with olivine sto- 
chiometry (65%) mixed with 1 and 20% carbon (abundances of 
50 and 15%, respectively). The abundance of amorphous carbon 
grains is of about 17%. Grains with sizes larger than 10 /^m con- 
tribute to 5% of the total emission in the 5-35 fim range. 

The crystalline phase: Emission features seen at 24, 28 and 
33 fim can be reproduced using Mg-rich olivine grains, with an 
abundance of 7%. The rather broad features around 28 fim are 
best matched with a additional contribution of enstatite (8%) 
and finally ~4% of yS-cristobalite are used in order to reproduce 
the shape around 20 fim. All these three crystalline species are 
mixed with 1 % carbon. 

Disk geometry: According to our best model, the dust is lo- 
cated at rather large distances from the star, with =53 AU and 
Tout =68 AU. The grain size distribution exponent is p = -3.6 
and we find a total dust mass of 1.2 x lO^^M© for grains with 
sizes smaller than lOOOjum. However, the results of the un- 
certainties estimation show that modeling the Irs spectrum of 
HD 169666 is an extremely degenerate problem. Given the low 
SNR over the entire spectral range, probability distributions for 
both Tin and rout are broad (cr ~20 and 15 AU, respectively). This 
means that the geometry of the debris disk cannot be firmly con- 
strained with existing data. 

Additional observations and comments: Our best model pre- 
dicts a flux at 70jum of 36mJy, which does not compare so 
well with the Mi ps measurement of 21.8+4mJy reported by 
iMoor et aH (l2009l) . The current best fit reproduces well the emis- 
sion features detected at longer wavelengths, however it does 
not match properly other, smaller emission features, especially 
those around 11 and 16jum. Among all the models tested by 
the genetic algorithm, we find another solution that reproduces 
these two specific features slightly better (but has a higher Xr)- 



For this model, the disk is located much closer to the star 
(Hn =13.4, rout=16.6AU versus 53 and ~68 AU), therefore the 
olivine grains are warmer and can better reproduce the emis- 
sion features around ll-16yum. To distinguish between both 
solutions, we refer to this one as the "warm" solution (with 
Tm-dx ~ 380 K) and to the original one as the "cold" solution 
(with Tmax ~ 215 K), and both fits to the data are shown in 
Figure[TT] As can be seen on the right panel of Fig.[TT] the warm 
solution does not provide a good fit to the 33 fim feature, which 
is the most prominent feature in the spectrum, and this model 
still over-predicts the Mips flux at 70 fim (32 versus 21.8 mJy). 
Compared to the cold model, the increases by ~7%. In the 
rest of the study, we will still consider the cold solution as the 
best fit model, since it provides the smallest and reproduces 
well the strongest emission features. But, this simple exercise 
underlines the crucial importance of SNR when modeling spec- 
troscopic observations. HD 169666 is the source with the lowest 
SNR once the stellar spectrum is subtracted, and this renders 
both the analysis and interpretations difficult. 

5.6. H D 98800 B 



Previous studies: As the title of I Verrier & Evans I d2008h sug- 
gests, HD 98800 B is an unusual debris disk. The system is com- 
posed of two bound spectroscopic bina ries A and B. Both sys - 
tems were estimated to be lOMyr old (ISoderblom et al.l 19981). 
With their Keck diffraction-limited observations. iPrato et al.l 
( 1200 ll) found that the dust belt surrounds the second spectro- 
scopic binary HD 98800 B, and that the dust belt is located be- 
tween ~ 2 and 5 AU. The structure of the dusty c omponent has 
also be en studied by several authors, for instance, iKoerner et al.l 
(I2OOOI) estimated the disk to be l ocated at Hn - 5 + 2 .5 AU with 
a width Tout - Hn of 13+8 AU. iFurlan et al] ( |2007|) also stud- 
ied the dust content around HD 98800 B with a two components 
scenario: an optically thin disk located between 1.5 and 2AU 
from the spectrosco pic binary and an optically thick disk at a 
distance of 5.9 AU. ILow et alJ (l2005h reported Mips fluxes of 
6.26+1.25 Jy and 2.1+0.42 Jy at 70 and 160 jum, respectively. 

The amorphous phase: A broad feature is seen at around 
10 fim, but no other emission features can be firmly detected 
in the Irs spectrum. Fe-bearing amorphous grains with pyrox- 
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Fig. 12. Fits to stellar subtracted spectra of HD 98800 B (fits in 
red, data in dashed black). Blue line shows the contribution of 
enstatite grains. 



ene stoichiometry accounts for 28% of the total abundance. For 
the olivine group, Fe-bearing grains represent 26% of the total 
abundance, and ~23% of Fe-free amorphous grains with olivine 
stoichiometry are used to reproduce the data. The abundance of 
amorphous carbon grains is 18%. Grains with sizes between 10 
and lOOOyum contribute to ~60% of the total emission in the Irs 
range. 

The crystalline phase: Even though no emission features at- 
tributed to crystalline grains are detected in the spectrum, the 
crystallinity fraction is not null. We tentatively detect Mg-rich 
olivine grains with a 0.2^^ 2'^ abundance. The presence of en- 
statite grains appears to be better constrained (5.2+ 1.5%). The 
contributions of both crystalline species peak around 19-20jum 
in order to try to reproduce the rather broad band between 18 
and ~21.5/im. Enstatite grains additionally contribute around 
28-30yL(m in order to match the other broad feature at the end 
of the spectrum. Figure[T2]shows the fit to the stellar subtracted 
spectrum as well as the contribution of enstatite grains to the 
total emission (in blue). We try to fit the spectrum without crys- 
talline dust species. The overall shape of the spectrum is repro- 
duced, however the best fit cannot account for the broad features 
around 20 and 30 ;um. Since the presence of crystalline grains is 
not supported by the detection of strong emission features above 
the continuum, we prefer to consider the presence of Mg-rich 
olivine and enstatite grains as tentative. 

Disk geometry: We find the dust to be located between 
Tin =1.3 AU and rout ~7AU, values consistent with previous 
analysis led on this object. The grain size distribution has an 
exponent p - -3, and we find a total dust mass of 6 x 10 ' Me, 
for grain with sizes smaller than 1000 /im. 

Additional observations and comments: Our model pre- 
dicts fluxes of 5.69 and 1.45 Jy at 70 and 160Atm, r espec- 
tively, to be compared with the values from iLow et aP (l2005l) 
of 6.26±1.25Jy and 2.1 +0.42 Jy at the same wavelengths. The 
small differences may be the consequence of even larger grains 
contributing to the emission, but obtaining a better match of the 
photometric measurement at 160yL(m is out of the scope of this 
study. 



5.7. [GBR2007JID8 

Emission in excess at mid-IR wavelengths for both 
[GBR2007]ID8 and [GBR200 7]ID9 (ID 8 a n d ID 9 here- 
after) was first identified by iGorlova et al.l (l2007l) while 
studying the NGC 2547 stellar cluster (d = 361 pc, 38+^-^Myr, 
lNavlor&.Teffriesll2006h . Both sources display departure from 
photosphe ric emission at ~4 a nd ~2p.m, respectively. Based 
on SEDs, iGorlova et al.l ( l2007h proposed three possibilities to 
explain the measured excess. First, disks are transiting between 
primordial toward debris disks, second a collision between two 
planetesimals occurred, producing large excesses and finally 
both objects could be background post-MS stars experiencing 
mass loss. Based on complementary UVES observations, ID 8 
was confirmed to be a member of the NGC2547 cluster, while 
ID 9 appears to be a background object rather than a member 
of the stellar cluster (Gorlova et al. in prep). In this study, we 
therefore focus only on ID 8. The Irs spectrum of this object 
shows a strong emission feature at 10 /im as well as at 16 and 
19/im. Unfortunately, as the source is faint, the Long Low 
module suffers from larger uncertainties, but a tentative feature 
at around ~ 24 jum is detected. 

The amorphous phase: Most of the dust in the disk around 
ID 8 is in its amorphous form (~98%). Fe-free amorphous grains 
with olivine stoichiometry clearly dominate the composition, 
with a total abundance of 61%. Fe-bearing amorphous grains 
from the pyroxene group are also required to reproduce the data, 
with an abundance of 26%. Amorphous carbon grains have an 
abundance of about 11%. Grains with sizes between 10 and 
lOOOyKm contribute to 4% of the total emission. 

The crystalline phase: The ~2% crystallinity fraction ob- 
tained for the dust composition is represented by Fe-rich and 
Mg-rich olivine grains (2 and 0.3 %, respectively), both mixed 
with 20% of carbon. The use of Fe-rich olivine grains instead 
of Mg-rich olivine mostly comes from the peak position of the 
feature around 24 //m. Given the low SNR in this spectral range, 
the detection of Fe-rich instead of Mg-rich olivine grains is con- 
sidered as tentative. The uncertainties show that Mg-rich olivine 
grains can also provide a satisfying fit to the data. 

Disk geometry: According to our best model, the dust is lo- 
cated close to the star with = 0.4 AU and rout = 0.5 AU. The 
grain size distribution is really steep, with a slope in -4.01 and 
the total dust mass is of about 2 x lO^^M® for grain with sizes 
smaller than 1000 /urn. The steep grain size distribution exponent 
can be explained by the rather unusual strength of the 10 //m 
emission feature, which requires that small grains dominate the 
emission. 

Additional comments: Interestingly, the spectrum of ID 8 
compares well with the spectrum of BD+20 307. Both spectra 
show a strong 10 //m emission feature, and a smaller flux level at 
longer wavelengths. Emission features from crystalline olivine 
grains are detected at 16, 19 and 24 jum. The main difference be- 
tween the two spectra is the emission feature at 1 1 jim on top 
of the emission feature from the amorphous grains for ID 8, and 
this is the reason why we obtain a good fit to the spectrum of this 
source and not for BD+20 307. 

6. Dust survival and origin 

In this Section, we discuss the possible origin of the dust belts 
using diagnostics for their luminosities and we investigate the 
main process that can be responsible for the evacuation of the 
smallest grains, in order to better constrain their survival and 
thus their origin. We also discuss the eventual presence of outer 
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Table 5. Observed (fobs) and maximal (/max) fractional luminosi- 
ties for our target list. 



Star /nax fobJfmnx Transient 



HD 11 3766 A 


7.2x 10- 


-2 


4.3x 10- 


-5 


> 10' 


/ 


HD 69830 


2.8x 10- 


-4 


5.9x10- 


-8 


> 10' 


/ 


BD+20 307 


5.2 X 10- 


_2 


2.6x 10- 


-7 


> 10' 


/ 


HD 15407 A 


7.0 X 10- 


-3 


7.0x 10- 


-6 


~ 10' 


/ 


HD 169666 


1.8x 10- 


-4 


1.6 X 10 


-4 


~1 


X 


HD 98800 B 


1.5x 10- 


-1 


2.8x 10- 


-4 


-700 


? 


ID 8 


3.2 X 10- 




4.0 X 10" 


-7 


> 10' 


/ 



dust belts located further away from the star in order to have a 
better description of the systems. 

6.1. Steady-state evolution or transient dust ? 

As explained in lWvatt et al.l (l2007h . the fractional luminosity of 
the debris disk /obs - Lir/L* is designed to be a proxy for the 
disk mass, in the optically thin case. If the observed debris disk is 
the result of a steady-state evolution, and assuming the disk has 
the same age as the host star, one can determine a maximum disk 
mass Mmax (and its corresponding maximal fractional luminosity 
/max) that can survive up to the age fare of the syste m. To compute 
/max, we start from Eqn. 20 of IWvatt et all (l2007h : 

/max = 0.58 X lO-'^r''\dr/r)D^:'Qf'e-'''M-J''Lf -'t;l,, (5) 

in which we inject the stellar parameters (fage, M* and L*), as 
well as our modeling results [r - r^n, dr = rout - Hn, both 
in AU) and use their prescribed values for the maximal size of 
planetesimals in cascade (Dc - 2000 m), planetesimal strength 
(2o = 200 J. kg"'), and planetesimal eccentricity (e = 0.05). The 
authors defined a threshold value of 100 for the ratio fohs/fmax- 
below this value, a steady-state evolution can explain the in- 
frared luminosity of the debris disk. Table|5]summarizes the val- 
ues we obtain for /obs, /max and their ratios. First, the case of 
HD 98800 B is difficuit to assess, the ratio b eing of about -700 . 
Given the relatively young age of the object. IWvatt et al] (l2007h 
concluded the collisional cascade has probably started recently, 
meaning the observed debris disk could possibly be the result of 
a steady-state evolution. Concerning HD 169666, we find a ratio 
of about ~1 for the best model ("cold" solution as discussed in 
Sect.l5.5|l and a ra tio of about 40 for the "warm" solution. As 
iMoor et all (l2009l) found a smaller value for (~ 4 AU), their 
ratio /obs //max is above the threshold value of 100 and the au- 
thors concluded the amount of dust is too important compared 
to what is expected for a steady-state model. Since the SNR is 
small for this source, we cannot easily conclude on the evolu- 
tionary status of this source since /max depends on the modeling 
results. Finally, the remaining five objects are too old for their 
high infrared luminosities to be explained by a steady-state evo- 
lution and are therefore considered as being in a transient phase. 
A dynamical event, responsible for the production of small dust 
grains, must have occurred so that we can observe such high 
fractional luminosities in the mid-IR. 

It is of interest to mention that most of the sources are either 
in binary systems or have known planets orbiting the central star 
The disks around HD 11 3766 A and HD 15407 A surround the 
primary of a binary system, and the dust belt around HD 69830 
is found to encompass the orbits of three Neptune-sized plan- 
ets. In the case of BD-i-20307, the dust surrounds the spectro- 
scopic binary, while the system of HD 98800 B is slightly more 
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Fig. 13. Radiation pressure efficiency (JSi-p) for grains around 
HD 15407 A, as a function of grain size s for amorphous 
MgFeSi04 (soHd black), MgFeSi206 (red dashed), MgSiOs 
(blue dotted), Mg2Si04 grains (green dotted-dashed), crystalline 
Mg-rich and Fe-rich olivine grains (solid magenta, and dashed 
cyan) and yS-cristobalite grains (dashed black). Shaded areas 
show regions where grains are gravitationaUy unbound to the 
star. 



complicated. The disk is located around the secondary of the 
HD 98800 A-B binary system and surrounds the spectroscopic 
binary HD 98800 Ba-Bb. The multiplicity of ID 8 remains un- 
known since the star was only recently discovered. Stellar com- 
panions or massive planets could be r esponsible for the stir - 
ring of smaller planetesimals (see e.g., iMatr anga et al 1 120101) . 
Therefore, we may be witnessing the consequences of dynam- 
ical instabilities triggered by massive objects in the vicinity of 
the central stars. 

6.2. Radiation pressure and time variabiiity 

Most of the debris disks of our sample are in transient phase, one 
may therefore wonder about the survival of the fresh, recently 
produced dust we observe. In Section|5] we have not discussed 
the effect of radiation pressure of the stellar radiation field on 
the small grains. This effect acts as an apparent reduction of the 
stellar mass and is the fastest way to evacuate //m-sized grains 
from dust belts. Using the absorption and scattering efficiencies, 
and assuming isotropic scatteri ng we compute the unitless effi- 
ciency ratio ySrp, as detailed in iBurns et all (Il979h for different 
grain sizes and different dust compositions. For a given grain 
size s and dust composition, ySi-p is the ratio F^p/Fg between ra- 
diation pressure and gravitational forces, respectively, and quan- 
tifies the competition between them. As both forces decrease as 
r"^, this ratio no longer depends on r. A p^, value greater than 
1 means that radiation pressure overcome gravitational forces 
and such grains will be evacuated from the system. In the range 
0.5 < y6,p < 1.0, grains are in the unbound regime with hyper- 
bolic orbits (see lKrivovll20I0h . and will be evacuated from the 
system as well. Figure[T3] shows an example of the /?,p -ratio for 
HD 15407 A, with s ranging between 0.1 and 10 //m for several 
dust species. An immediate result is that grains smaller than typi- 
cally 3.5^ jjm are rapidly ejected from the system {J3i-p > 0.5). If 
a grain is subjected to a fS^p ratio greater than 0.5, then the eccen- 
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tricity of its orbit will increase and the grain will quickly be evac- 
uated from the system (comparable to the dynamical timescale). 
In the case of HD 15407 A, the evacuation of grains smaller than 
~ 3.5 fim is expected to be very efficient, of the order of a few 
months. 

Among the sources studied in our sample, grains around 
HD 98800 B and ID 8 are not su bjected to evacua t ion by ra- 
diation pressure. As discussed in iBeichman et all (1201 ll) for 
HD 69830, only grains with sizes in a narrow range (0.4 < s < 
0.7 yum) are affected by radiation pressure. Concerning the spec- 
troscopic binary BD+20 307 (M* = 2.2M0), grains with sizes in 
the range 0. 3-0.7 ;um (0.4-0.9 //m, depending on the composi- 
tion) are expected to be evacuated from the system. This means 
that grains smaller than ~0.3 fim may survive to radiation pres- 
sure around the binary system, and these small grains can con- 
tribute to the emission of the strong and sharp lOyum feature. 
Overall for both HD 69830 and BD+20 307, since the ranges of 
grain sizes with ySip > 0.5 are narrow, we consider this effect to 
be of second order for the purpose of this study. Additional con- 
straints are necessary to lower the degeneracies in the problem 
in order to refine the model and assess the influence of radia- 
tion pressure. For the four sources discussed here, the survival 
of the dust belt in time (and hence their detection in the mid-IR) 
is less constraining on the mechanisms responsible for the pro- 
duction of small yum-sized grains, compared to the other three 
sources. Since the timescale for evacuation by P-R drag is much 
longer (fp_R ~ 400y6|:p' n M^' yrs, with r in AU and in Mq, 
IWvattll2008l) than the timescale for evacuation by radiation pres- 
sure (< 1 yr), these small grains can survive much longer being 
gravitationally bound to the central star. 

In the following, we discuss the effect of radiation pressure 
for dust grains in the debris disks around the remaining three 
sources (HD 15407 A, HD 113766 A and HD 169666). For each 
sources, we discuss which grains are in the unbound regime, de- 
pending on their sizes s. It is however out of the scope of this 
paper to include the effect of radiation pressure in the modeling 
approach. It would indeed require to have robust prescriptions 
about complex, and time-dependent, dynamical processes such 
as dust evacuation and dust fragmentation. Instead we discuss 
timescales about variability and subsequent constraints on the 
production of small //m-sized dust grains. 

6.2.1. Time variability of HD 15407 A 

As shown in Fig.[T3] amorphous grains smaller than ~ 3.5^;um 
are short-lived around HD 15407 A. Interestingly, we find that 
some small y6-cristobalite grains, mixed with 1 % carbon and with 
sizes in the range 0. 1-0.25 //m, can survive to radiation pressure. 
Larger grains or grains mixed with higher carbon fractions have 
ySip greater than 0.5. According to our results (Sect.|53) silica 
grains with an addition of 1% carbon provide the best match 
to the data. This finding slightly relaxes the constraints on the 
production of silica grains. However, we could not obtain a sat- 
isfying fit to the data (especially to the emission features from 
jS-cristobalite) without amorphous grains smaller than 3 fim. 

To assess the time variability of the debris disk around 
the primary, we gather photometric observations that were per- 
formed in the last decades for this object. Table|6] summarizes 
our findings and reports the instrument used and wavelengths 
of observations as well as measured fluxes. The column labeled 
"PiRs" reports flux predictions based on the Irs spectrum and 
their uncertainties. To obtain these predictions and their uncer- 
tainties, we convolve the Irs fluxes to the filters widths of other 
instruments. The photometric observations cover a range of sev- 



eral decades. Unfortunately, in some cases determining an ac- 
curate epoch for the measurement is not trivial. All of the rele- 
vant facilities and instruments listed in Tab.|6]observed the vicin- 
ity of HD 15407 A several times during their mission and the 
quoted flux densities in the specific catalogues have been de- 
rived as a combination of the photometry extracted from indi- 
vidual scans or observations. However, we know that the MSX 
Galactic Plane survey observations were obtained between May 
1996 and January 1997, while AKARI was active between May 
2006 and August 2007. The IRAS mission was accomplished in 
1983. Utilizing the SCANPI software we find that the relevant 
IRAS scans were performed between the 17''' of August and the 
15'''of September 1983. In the case of WISE, the observations 
were performed the 8''' and 9''' of February 2010. 

The measured fluxes over almost 3 decades compare very 
well with the fluxes derived from the Irs spectrum (observed on 
the 9''' of October 2008), especially at 9 fim where the emission 
feature from small y6-cristobalite grains peaks (960+57 versus 
927+47 mJy). At least one year passed between the AKARI and 
Spitzer observations, a timescale to be compared with a radia- 
tion pressure timescale of a couple months. We therefore con- 
clude that the system has been collisionally active over several 
decades and that the population of grains smaller than 3 yum is 
being efficiently replenished over time. 

6.2.2. Time variability of HD 1 13766 A 

In the disk around HD 11 3766 A, all grains (amorphous and 
crystalline) smaller than ~3.5;um should in principle be evac- 
uated from the system on a timescale of months. This result 
contrasts with the strong and sharp emission features detected 
in the spectra, arising from amorphous grains as probed by 
the 10/im feature and from crystalline olivine grains (e.g., at 
llyum). As for HD 15407 A we gather photometric observa- 
tions of HD 1 13766 A, and compare them with the Irs spectrum 
convolved with the corresponding photometric filters. Table|6] 
summarizes the measurements and predictions. Again, obtain- 
ing dates of observations is difficult, however the IRAS mission 
was accomplished in 1983, AKARI was active between May 
2006 and August 2007, and Mips observations were performed 
the 23"' of February 2004. WISE observations were obtained in 
January (25-26"') and July (23-24"') 2010. The photometric ob- 
servations compare very well with the convolved fluxes, mean- 
ing the emission in excess around HD II 3766 A should have 
been stable for the last 25 ye ars. In addition to t heir VLTI/Midi 
interferometric observations, ISmith et aP (l20I2l) also presented 
VLT/VisiR spectroscopic observations of HD 11 3766 A, in the 
8-13yL(m spectral range. Interestingly, the shapes of both Visir 
and Irs spectra compare very well with each other (even though 
the new observations return a flux level that is of about 1.3 times 
smaller compared to the 2004 Spitzer observations). The pres- 
ence of emission features, and hence of warm //m-sized grains 
indicates that these dust grains must be replenished quickly, on a 
timescale of a few orbits. This strongly points towards on-going 
collisions in the disk. 

6.2.3. Time variability of HD 1 69666 

Computing the efficiency ratio /3,p for grains of different com- 
positions (MgFeSi04, MgFeSi206, MgSiOa and Mg2Si04) we 
find that grains smaller than 4jjm are expected to be short- 
lived in the disk around HD 169666. And yet emission fea- 
tures arising from small crystalline olivine grains are detected 
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Table 6. Photometric observations (columns "Flux") of HD 15407 A, HD 1 13766 A and HD 169666 and flux predictions from their 
Irs spectra convolved with the corresponding photometric filters (columns "Pirs")- Second line of the Tabel give the epochs of the 
Irs observations (day /month/year). 









HD 15407 A 






HD 11 3766 A 




HD 169666 










(09/10/2008) 






(01/03/2004) 




(02/07/2005) 




Facilities 


A 


Epoch 


Flux 


Pirs 


Epoch 


Flux 


Pirs 


Epoch 


Flux 


Pirs 




[jjm] 




[miy] 


[mIy] 




[mJy] 


[mJy] 




[mJy] 


[mJy] 


MSX 


8.28 


1996-1997 


808 ±49 


894 ±45 














AKARl/lRc 


9.00 


2006-2007 


960 ±57 


927 ±47 


2006-2007 


1308 ±92 


1439 ±72 


2006-2007 


344 ± 18 


333 ± 17 


WISE 


11.56 


Feb. 2010 


641 ±30 


673 ± 34 


Jan.-M. 2010 


1161 ±53 


1372 ±71 


Feb.-Mar. 2010 


172 ±8 


182±9 


IRAS 


12.00 


1983 


1050 ±63 


891 ±45 


1983 


1590 ±79 


1746 ±89 


1983 


232 ± 16 


242 ± 12 


AKARVIrc 


18.00 


2006-2007 


497 ± 33 


445 ± 22 


2006-2007 


1428 ± 87 


1504 ±79 


2006-2007 


130 ±23 


107 ±6 


WISE 


22.09 


Feb. 2010 


400 ± 24 


384 ±19 


Jan.-M. 2010 


1656 ±96 


1554 ±81 


Feb.-Mar. 2010 


103 ±6 


99±5 


Spitzer/Mips 


23.67 








Feb. 2004 


1459 ±30 


1431 ±74 


Jun. 2005 


90±4 


94 ±5 


IRAS 


25.00 


1983 


432 ±39 


402 ± 20 


1983 


1760 ± 123 


1664 ±87 


1983 


101 ± 13 


112±6 



at longer wavelengths. iMoor et al.l (|2009|) already discussed the 
time variability of this debris disk. First, they compared two 
Irs spectra, obtained on the 2'"' of July 2005 and on the 6''' of 
June 2008 (we model the first one in this study as it includes 
the 2" '' order of module SL, between 5.2-8.7 jum). iMoor et aP 
(l2009l) concluded that, in a 3 year span, no significant differ- 
ences are seen between both spectra, and emission features can 
still be detected at longer wavelengths. In a similar way as for 
HD 15407 A and HD 113766 A, we gather photometric observa- 
tions of HD 169666 and compare them with the Irs spectrum 
convolved with the photometric filters. Last two columns of 
Table|6]show both observed and predicted fluxes. The strongest 
constraint on the time variability comes from the IRAS obser- 
vations of 1983, for which the observed and predicted fluxes 
compare extremely well ( Fv.obs=232 mJy against Pirs =242 mJy 
at 12 yum, and Fv,obs=101 mJy against Pirs=1 12mJy at 25 /im). 
Spitzer/Mips observations were obtained on the 23''' of June 
2005, and WISE photometric observations were performed be- 
tween the 21" of February and 6''' of March 2010. As for the 
previous two sources, there must be an on-going process in the 
disk that continuously replenishes the small dust grain popula- 
tion, as the transient phase has been stable for at last 25 years. 

6.3. Presence of outer, cold belts 

The dust detected around the debris disks we study is transient, 
except maybe in the cases of HD 98800 B (given its relatively 
young age) and HD 169666 (given the uncertainties on the lo- 
cation of the dust belt because of the low SNR in its spectrum). 
Two main scenarios have been invoked in the literature to explain 
the origin of the transient dust. First, the dust may be produced 
directly close to the star This scenario can be the consequence of 
several events. Transient dust may originate from a massive col- 
lision between two large bodies being disrupted close to the star. 
Even though the possibility of a steady-state asteroid belt is most 
likely ruled out by the high fractional luminosities of the disks, 
transient dust may be produced in an asteroid belt that has been 
recently excited and where a collisional cascade started. The sec- 
ond main scenario involves the presence of an outer belt dynam- 
ically unstable feeding the inner regions of the disk. Bodies from 
the outer belt are being scattered inward and rele ase small grains 
as they collide or sublimate. As mentioned in iBeichman et al.l 
(l2005h for HD 69830, a "super-comet" from the outermost re- 
gions may have been captured on a low eccentricity orbit in the 
inner regions and is being evaporated and broken into smaller 



bodies. However, detecting an outer belt further away from the 
star does not enable us to disentangle between the two scenarios, 
since the outer belt may be stable in time and not responsible for 
the production of the transient dust. 

The most straightforward way to infer the presence of cold, 
outer dust belts is via far-IR, or mm observations. For one source 
in our sample (HD 11 3766 A), we can already mention that an- 
other dust belt is likely to be located further away from the star 
Spitzer/Mips observations at 70 yum reveal an excess of ~ 1 20 mJy 
compared to the predicted flux of our best model (270 mJy). 
As Irs observations are unresolved, we cannot quantify by how 
much the outer belt contribute to the flux measured by the spec- 
trograph. In this study, we assume the dust is in a single, contin- 
uous disk. Consequently, the flux contribution of the outer belt in 
the Irs range translates into rather large values for rout (12 AU). 
This can be explained as a need for a contribution of colder 
grains. The extension of disks, in the formalism of our model, 
may be a possible diagnostic for the presence of outer dust belts. 
According to our findings, HD 113766 A is not the only source 
for which we find an extended disk. The width of the disk around 
HD 15407 A is of about -19 AU, which also suggests the need 
for cold dust grains in order to reproduce the unresolved spec- 
troscopic data. Further observations are required in order to dis- 
tinguish between the two possibilities of either an extended disk 
or two spatially separated dust belts. Even though the Mips ob- 
servations at 70 ijm do not reveal emission in excess compared 
to our model, we discussed the possibility of a two belts sce- 
nario around BD+20307 (Sect. l5.3l l. The justification being that 
the disk is rather extended (18 AU) and that we observe a gradi- 
ent in the dust composition: the warm inner regions probed by 
the lOjjm emission feature seem to contain pristine dust, while 
emission features from crystalline olivine grains are detected at 
longer wavelengths and may arise from another dust belt. With 
the assumption of one continuous, homogenous disk and even 
allowing olivine grains to have low temperatures (mixed with 1 
or 5% of carbon), we cannot achieve a satisfying fit to the 10 //m 
feature. 

The four remaining objects in our sample (HD 69830, 
HD 169666, HD 98800 B, and ID 8) could be modeled success- 
fully with the assumption of a single dust belt. For the first three 
objects complementary observations at longer wavelengths do 
not show any emission in excess, except for the 160/^m Mips 
photometric point for HD 98800 B. In this case, including larger 
grains in the model may compensate for the small difference. 
The absence of far-IR emission in excess for these sources points 
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towards the scenario of a dust production mechanism directly 
close to the star. However, determining which declinations of 
this scenario (massive collision, grinding of an asteroid belt, 
super-comet) is the correct one is difficult. In the following 
Section, we will investigate if the results from spectral decom- 
position can provide additional information on the origin of the 
dust. 



7. The origin of crystalline grains 

To have a better understanding of these transient debris disks, 
one can investigate the origin of crystalline grains. There are 
several possibilities to crystallize amorphous grains: (i) thermal 
annealing, where grains are continuously subjected to temper- 
atures high enough (or timescales long enough at intermediate 
temperatures), for their internal structure to be re-an^anged, (ii) 
shock-induced melting, where grains are subject to high temper- 
atures for short amounts of time, due to collisions. Finally, (Hi) 
crystalline grains detected in the disks may originate from the 
parent bodies. The parents bodies are either asteroids progres- 
sively grind down in a belt, cometary bodies that were scattered 
inwards in a LHB-like event, or large planetesimals that collided 
and were completely disrupted. In the latter case, either colli- 
sions produce enough heat for crystallization, or crystals were 
already formed through differentiation in parent bodies and are 
released through collisions. Scenarios (ii) and (Hi) are intimately 
connected since they both involve collisions between planetes- 
imals or smaller bodies, but in the first case, collisions release 
enough energy to melt the fragments, while in the second case, 
the fragments are not altered by the collisions. As debris disks 
are gas-poor, we do not considered "present-day" gas-phase con- 
densation as a possible significant source of crystalline grains. In 
the following paragraphs, we present and discuss which diagnos- 
tics on crystallization can help us better understand the origin of 
these grains. 



7.1. Where's the enstatite ? 

The spectra in our target list do not display any strong, recog- 
nizable emission features that can be attributed to enstatite crys- 
talline grains (e.g., at 9.3/im). This translates into overall small 
abundances for this particular dust species, while enstatite has 
been commonly observed around a large variety of primordial 
disks around Class II object s, with different ages, stellar masses 
and star forming regions ( Bouwman et al. 20081 lJuhasz et al.l 
120101 lOlofsson et al.ll201flHSargent et al.ll2009bl) . One result of 
these studies being that enstatite seem to dominate the crys- 
tallinity fraction of the inner regions (~ 1 AU for TTauri stars), 
while forsterite dominates the outer regions of the disks (< 10- 
15 AU for TTauri stars). Concerning debris disks, so far only 
one object seems to harbor signifi cantly more enstatite grains 
than forsterite grains (HP 165014. iFuiiwara et aDl2010h . This 
contrast between gas-rich, massive disks and gas-poor debris 
disk is intriguing. For the objects we study, if the observed dust 
was the product of a collisional cascade of planetesimals formed 
within the first AU of the star, we may expect to detect more 



enstatite dust grains. This is one possible scenario mentioned by 
iFujiwara et al.l (l2010l) . who speculate that the debris disk around 
HD 165014 may be the aftermath of a collision between small 
bodies and the surfac e of a Merc ury-like planet (which shows 
evidence for enstatite. ISprague & R oush 1998). The overall low 
abundances of enstatite, and the lack of strong emission features, 
pratically rule out such events for the sources in our sample. 



The observed crystalline grains may therefore originate from 
bodies that were formed at larger distances from the star (where 
enstatite is scarcer), or the observed crystalline grains were re- 
cently formed inside the disks. Both possibilities fit rather well 
the picture of low enstatite abundances. If the debris disks mod- 
eled in this study are the results of collisions between bodies 
originally formed further away from the star, the lack of en- 
statite grains compares well with observations of younger disks. 
On the other hand, if the observed crystalline grains are formed 
at present days in the debris disk, laboratory experiments can 
help us understanding the low abundances of enstatite. For in- 
stance, wit h their annealing experi ments of amorphous grains 
at 1000 K, [Thompson etM ( l2002l) have shown that enstatite 
is difficult to form at these temperatures, and forsterite grains 
are preferentially formed (as probed by X-ray powder diffrac- 
tion). This result echoes the respective activati on energies for 
annealing of enstatite and forsterite grains (see lGaiill2010l and 
references therein). Several laboratory experiments have shown 
that activation energies for annealing are larger for enstatit e 
grains compared to forsterite grains (e.g., iFabianet al.ll2000l) . 
Consequently, if present day thermal annealing is the main pro- 
duction source of crystalline grains in debris disks, we can ex- 
pect to have higher abundances of crystalline olivine grains com- 
pared to enstatite. Interestingly, our finding of low abundances of 
enstatite echoes the study of crystallization processes during the 
violent outburst around EXLupi. Abraham et al . (2009) found 
evidences for production of forsterite, but not of enstatite, within 
the first AU of the disk, a result later on confirmed by VLTI/Midi 
observation (lJuhasz et al.ll201 ll) . 



7.2. The iron content 

In order to distinguish between the different scenarios on the ori- 
gin of crystalline grains, the petrology of the olivine dust grains, 
and more specifically the iron content of these grains, can pro- 
vide new insights on the crystallization processes at stake in 
these disks. This means, we have to be confident in our results 
about the Fe content of the olivine grains. This study relies on the 
laboratory measurements for different olivine samples. We in- 
sure to use data with as little contamination as possible (aerosol 
data). Additional effects may come into play, such as the shape 
of the grains, that can possibly shift peak positions (iKoike et al.l 
1201 Oh . but also po rosity or crystal line grains being embedded 
in other materials (iMin et al.ll2008h . Concerning the first point, 
many olivine grains are injected in the aerosol apparatus, with 
various sizes and shapes. Plus, we do not make use of any scat- 
tering theory, such as the DHS the ory that may oversimplify 
some aspects of the problem (Mutsch ke et al.l2009h . Concerning 
the last two additional effects, we unfortunately cannot take them 
into account as this will result in too many degeneracies for an 
already complex modeling. We would like to stress again the 
need for high angular resolution observations to constrain some 
crucial parameters in the modeling. To conclude, we are con- 
fident about the detection of Fe-iich olivine grains in several 
debris disks (see Sect. l4.2l i. given the available data and obser- 
vations and the assumptions made on the disk geometry (e.g., 
continuous dust belt). 

This finding of Fe-rich olivine grains contrasts with observa- 
tions of Class II disks, for which most studies concluded on the 
absence of Fe-rich crystalline silicates. In their sam ple of Herbig 
Ae/Be stars, based on the 9.4yL(m emission feature. lJuhasz et aH 
( l20IOh reported a possible fraction of 10% of Fe in crystalline 
pyroxene grains. However, the authors do not report any shift to 
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longer wavelengths for emission features associated with crys- 
talline olivine grains (e.g., 19 or 24 yum), which indicates an over- 
all small fraction of Fe in olivine grains. To shed light on these 
petrological differences of olivine grains, we discuss the results 
of annealing experiments in laboratory. If amorphous dust grains 
are exposed to high temperatures for a significant amount of 
time, internal re-arrangement processes are activated in the ma- 
terial and a periodic order appears in the originally disordered 
network (see lHenning|20 1 Ol) . Crystallization is a process that de- 
pends on two main environmental parameters: temperature and 
time. Amorphous grains can crystallize either at intermediate 
temperatures if the exposure time is long enough, or on short 
timescales if the temperature is high enough. In the following, 
we summarize some of the relevant results from previous studies 
about thermal annealing, with respect to the influence of these 
tw o parameters on dust petr ology. 

iNuth & JohnsonI (l2006h investigated the production of Fe- 
bearing crystalline grains via thermal annealing. According to 
their results, both end-members of magnesium and iron crys- 
talline grains (forsterite and fayalite for the olivine group, re- 
spectively) can be formed via thermal annealing, if the temper- 
atures are high enough (~ HOOK and 1400 K, for magnesium 
and iron, respectively). However, as the temperatures required 
to form crystalline iron silicates are higher than temperatures re- 
quired to form crystalline magnesium silicates, the authors argue 
the crystalline iron grains will evaporate on timescales of several 
m onths, while magnesiurn grains can survive longer According 
to lNuth & JohnsonI ( l2006l) . this timescale issue is the reason why 
iron crystalline grains are not detected in young protoplanetary 
disks: the grains are destroyed before they can be transported 
outwards and survive. Since iron end-member crystalline grains 
(e.g., fayalite) were not required to fit our data, this tempera- 
ture threshold (1 100 vs. 1400 K) is not directly applicable to our 
problem. Still, iNuth & JohnsonI (l2006b discuss thermal anneal- 
ing of Fe- and Mg-silicate aggregates and suggest as a first ap- 
proximation that annealing conditions for such aggregates might 
fall in the 11 00-1 400 K temperature range. Therefore tempera- 
tures of 1200-1300K might be sufficient to produce the Fe-rich 
(Fe / [Mg + Fe] ~ 20%) crystalline grains that we detect. One 
should still keep in mind that the exact conditions for anneal- 
ing of such mixed aggregates are yet to be better studied, which 
ren ders the discussion c hallenging. 

iMurata et al.l (l2009l) have conducted annealing experiments 
of Fe-bearing amorphous silicates during 600 min at tempera- 
tures of about ~1000K, at a pressure of one atmosphere. The 
oxygen partial pressure in the furnace was controlled, to avoid 
possible oxidation of Fe^^ cation. They conclude that the Mg 
fraction (=Mg/[Mg+ Fe]) of the crystalline phase decreases 
with time, as the amorphous phase is annealed. At the end of 
the experiment, the crystallinity fraction in the original material 
reached ~ 80% and the Mg fraction was close to 0.8. Such frac- 
tion is quite similar to the Fe-rich olivine grains we use in this 
study. We can speculate that annealing over long period of time 
will lower the Mg fraction to match the one of the inital amor- 
phous material. 

7.3. Collisions and shock-induced melting 

Transient dust observed around debris disks may be the con- 
sequence of collisions, where either km-sized bodies are con- 
tinuously being grind down, or larger planetesimals are com- 
pletely disrupted after a collision at high relative velocity (vrei ~ 
lOkm.s '). One can wonder how much energy these disruptive, 
or grinding collisions release, and by how much the temperature 



of the resulting fragments can increase. Using laboratory shock 
experiments as well as numerical simulations iKeil et all (Il997h 
studied the efficiency of shock heating by collisions (with Vrei 
ranging from 3 to 7 km.s"') for planetesimals with sizes between 
10 and 1000 km. Their conclusions were that impacts do not sig- 
nificantly heat up the fragments, mostly because of the low sur- 
face gravity of the planetesimals. For instance, they find a glob- 
ally averaged increase of about 50 K for a collision involving 
a 1000 km-siz ed body. But the planetesimals considered by the 
simulations of iKeil et alJ (Il997h were non-porous bodies. Since 
many ast eroids in the Sol ar System have non negligible porosity 
fraction dBritt et al.ll2002l) . more refined models should consider 
the ca se of porous bodies. As demonstrated by iDavison et al.l 
(1201 Ol) . including porosity in simulations can make a significant 
difference. The main conclusions from their study were that be- 
fore disruption, shock-induced heating is relatively minor and 
localized, but if the collision is disruptive then heating can be 
significant, with more than ~ 10% in mass of the parent body 
is melted. Collisions between porous planetesimals of compa- 
rable sizes can melt almost the entire mass of the two bodies, 
even at collision velocities smaller than 7km.s '. Since many 
uncertainties remain, applying the results of these simulations to 
our results is difficult. Nonetheless, we cannot rule out disruptive 
collisions as a source of newly formed crystalline grains. As the 
maximum temperature reached during such events is unknown, 
we can hardly address the question of the petrology of oUvine 
grains when produced by colUsions. 



7.4. Crystallization in differentiated bodies 

Interiors of non-porous planetesimals are known to be regions 
of high temperatures. The formation of planetary cores requires 
the melting of accreted bodies, and consequ ent differentiation. 
As discussed bv lMoskovitz & Gaidosl(l201 ll) . the heating source 
can be the decay of short-lived radionuclides ^^Al and ^"Fe, 
which can result in temperatures high enough for the melting 
to take place. Differentiated bodies will then be the result of 
gravitational segregatation of molten silicates and metals within 
the plan etesimal. According to studies of iron meteorites (e.g., 
Schersten et al 1 120061 lOin et"ani2008i) . such meteorites are the 
consequence of accretion, differentiation as well as crystalliza- 
tion of their original parent bodies, suggesting that such event 
most likely happened in the early Solar System, when the decay 
of ^''Al was significant enough to contribute to the heating. Even 
though the question of Fe-enrichment of olivine grains in diff'er- 
entiated bodies is at present day a poorly constrained question, 
we hypothesize that Fe-rich olivine grains may be produced in- 
side young differentiated bodies. Disruptive collisions between 
these bodies can then release the crystalline grains ( and not al- 
tering them in non-porous collisions as studied by iKeil et all 
1 19971) . an hypothesis that could account for the origin of the 
thermally metamorphosed LL chondrites, that contain Fe-rich 
olivine grains dNakamura et al.ii201 li) . 



7.5. Individual source discussions 

Here, we discuss for each source which scenario can best explain 
the inferred dust mineralogy, based on the several diagnostics 
presented previously. 
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7.5.1. HD 113766 A 

According to our best model, is located close to the sub- 
limation radius (Td.max ~ 1200 K). Fe-rich olivine grains are 
found to be more abundant than Mg-rich olivine grains. Since 
olivine grains produced in Class II disk are preferentially Mg- 
rich, we are not observing the mineralogy of parent bodies in- 
herited from the primordial disk (e.g., comets, see Sect. 18. Il l as 
the origin of the dust around HD 11 3766 A. Additionally, from 
IRAS observations, we know that the transient phase has last 
for at least 25 years, suggesting a continuous replenishment 
of yum-sized grains. Since disruptive collisions of bodies inher- 
ited from the primordial disk do not support our findings, we 
first propose that we may be witnessing the continuous grinding 
of several differentiated bodies (see Sect. 17. 41 and 18.3b . Fe-rich 
olivine grains could have been produced inside large planetes- 
imals and only then released in the exo-zodiacal belt as a con- 
sequence of disruptive collisions. Second, given the relatively 
high dust temperatures we find for our best model, one can 
also hypothesize that crystalline grains are pr oduced via ther- 
mal an neaUng, in a similar way as discussed bv lNuth & Johnsoiil 
(l2006l) . Crystalline Fe-rich silicates may be formed at temper- 
atures in the range 1100-1400K from Fe-, Mg-silicate aggre- 
gates. The eventual timescale issue for the survival of these 
grains is then comparable to the timescale of evacuation by 
radiation pressure for the (sub-) //m-sized olivine grains. With 
this picture in mind, we do not have to assume that all the 
observed yum-sized crystalline grains are continuously replen- 
ished from the parent bodies. Instead, we can suggest that a 
fraction of the small amorphous grains is constantly being an- 
nealed into Fe-rich olivine grains. The tentative detection of en- 
statite (0.8+0.8%) fits in both scena rios. Very few Solar System 
chondrites contain enstatite (~ 2%. iNorton & Chitwoodll2008h . 
which points toward comparable composition of the parent bod- 
ies around HD 1 13766 A. In the case of present-day crystalliza- 
tion, the timescale for production of enstatite may be longer than 
the timescale for evacuation of the grains. 

7.5.2. HD 69830 

The dust belt is located slightly further away from the 
star compared to the sublimation radius of silicate grains 
(T±max ~ 700 K) . There fore, thermal annealing as described by 
iNuth & Johnso rJ (I2006h cannot account for the abundances of 
Fe-rich olivines we find. However, radiation pressure has very 
limited efifect on sub- yum-sized grains, meaning t he grains sur- 
vive m uch longer in the disk. As explained by iMurata et al.l 
(l2009l) . and since survival timescales are no longer an issue, 
long exposure times to intermediate temperatures can produce 
Fe-rich crystalline grains, annealed from Fe-bearing amorphous 
grains that dominate the amorphous phase according to our re- 
sults. We could possibly be witnessing the effect of Mg/Fe frac- 
tionation in olivine grains. However, one open question would 
be to compare the timescales for fractionation and collisions. 
According to our result the s ystem can be de scribed as a col- 
lisional case ade (p = -3.54. lDohnanvilll969l) . suggestmg nu- 
merous collisions at all grain sizes. In a debris disk, relative ve- 
locities are too important for coagulation between grains to take 
place, and collisions most likely result in fragmentation of the 
colliding bodies. Then, how would fragmentation affects the ef- 
ficiency of fractionation? One may hypothesize that the energy 
released in the collision may accelerate the fractionation process 
by increasing the grains temperatures, but we cannot conclude 
on such conjectures without further investigations on that mat- 



ter On the b a sis of the activation energy for Mg2Si04 smoke in 
iFabian et al.l (120001) . the characteristic crystallization timescale 
at 700K is of about 3000 yr However, at 600K the timescale in- 
creases up to 30 Myr, to produce pure forsterite. Since the activa- 
tion energy of iron-rich silicates should be higher than Mg-rich 
silicates this timescale can increase significantly for iron-bearing 
silicates. This renders the possibility of thermal annealing as the 
main production source for crystalline grains quite unlikely. 

As suggested by the absence of an outer belt, the origin 
of the transient dust around HD 69830 is most likely the re- 
sult of collisions close to the star Since Fe-rich olivine grains 
are unlikely to be formed in Class II disks, we can potentially 
rule out that the parent bodies are inherited from the primordial 
disk. Therefore the possibility of an asteroid belt being slowly 
ground down may not fit very well in this picture. Comets like 
Hale-Bopp or Wild 2 are thought to have formed in the outer 
regions of the Solar System, and the crystalline olivine grains 
observed in both comets are Mg-rich (see Sect. 18. U . In the sce- 
nario of a super-comet be i ng br oken in smaller bodies, as de- 
scribed in lBeichman et al.l (l2005l) . we would therefore expect to 
preferentially detect Mg-rich olivine grains, formed in the pri- 
mordial disk and transported outside afterwards. The petrology 
of the olivine grains we infer in our best model does not sup- 
port this scenario. Shock-induced melting could contribute as 
a source of crystalline grains, but it would imply that crystal- 
lization via melting strongly favors the production of Fe-rich 
crystalline grains, a question that is poorly constrained at the 
moment. However, we can hypothesize that the transient dust 
originate from differentiat ed bodies, inside which olivines grains 
we re enriched in Fe (e.g.. lNakamura et al.ll20lTI) . As mentioned 
by iKeil et all jI997l) . if such planetesimals are non-porous (a 
reasonable assumption for differentiation to take place), their 
breakup into smaller fragments should not alter the mineralogy 
too much since shock-induced melting then becomes negligible. 
Interestingly, the system is the host of three Neptune-sized plan- 
ets. The observed dust belt may be the result of a dynamical 
instability, possibly triggered by planet migration for instance. 
Such an instability may have resulted in the collision and dis- 
ruption of large planetesimals. 

The small abundance of enstatite may correlate with the 
abundance of siHca (Mg2Si04 + SiOz ^ 2 MgSiOs). Production 
of enstatite via thermal annealing may not fit well in this picture, 
since the timescale for annealing at these intermediate temper- 
atures (700 K at maximum) is most likely longer than the colli- 
sional timescale. The finding of low abundance of enstatite com- 
pares well with the overall scarcity of enstatite chondrites in the 
Solar System. 

7.5.3. IDS 

The origin of crystalline grains around ID 8 is more diflicult to 
discuss because of the low SNR at longer wavelengths. This 
renders the spectral decomposition more uncertain on the ex- 
act crystallinity fraction as well as on the iron fraction in olivine 
grains. The uncertainties from the MCMC run show that the de- 
tection of Fe-rich over Mg-rich olivine grains is less convinc- 
ing compared to the previous two sources. We therefore prefer 
not to conclude on the exact petrology of the crystalline olivine 
grains in this system. Still, one of the interesting characteristic 
of the modeling results is the very steep grain size distribution 
(p = -4.01, for s up to 1000 /im). This suggests the system is at a 
peak of collisional activity, therefore producing large amounts of 
sub-jum-sized grains, a majority of them being amorphous. It in- 
dicates that whatever the parent bodies are (asteroids, comets or 
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planetesimals), they should have had an overall low crystallinity 
fraction. 

7.5.4. HD 15407 A 

The spectrum of HD 15407 A does not show emission features 
associated either with olivine or enstatite grains. But the spec- 
trum shows several emission features associated either with 
cristobalite or coesite (see Sect. l5.4l i. The exact nature of the 
Si02 polymorph is difficult to constrain since they both show 
an emission feature around 16/^m. Additionally, it is possible 
that other polymorphs (e.g., non-crystalline silica) are present 
in the debris disk around HD 15407 A, but inferring their pres- 
ence is a degenerate problem given the strong similarities be- 
tween the emissivity profiles, especially for the strongest emis- 
sion feature, around 9 yum. The /j-cristobal ite form of silica w e 
use is close to the annealed silica studied bv lFabian etan(l2000h . 
which is the result of an annealing experiment, for 5 h at a tem- 
perature of 1220 K. This indicates that the grains have experi- 
enced r ather high temperatures, but as discussed in lSargent et al.l 
(I2009ah . once formed cristobalite grains must be cooled quickly 
enough to retain their crystalline structure. If not, the grains are 
expected to revert to the lower temperature polymorphs of silica 
(j6- and a-quartz). Formation of coesite, instead of j6-cristobalite, 
is not a straightforward process either since it requires both high 
temperatures and pressure^ Such conditions may be fuUfilled 
in a collisional environment. It is noteworthy to remark that in 
the case where pressures are not high enough to form coesite, 
y6-cristobalite should be formed preferentially. 

The mineralogy of the dust around this source is puzzling. 
No crystalline silicate grains are detected (olivine or enstatite), 
but a significant fraction of SiOa is required to model the spec- 
trum. According to our best model, the maximum temperature is 
of about 1530 K, a temperature high enough to anneal silicates 
grains and produce for example crystalline olivine grains. Even 
though small grains are short-lived around HD 15407 A because 
of radiation pressure, if the main source of silica was thermal 
annealing, we should expect to detect features associated with 
olivine grains (as in the case of HD 1 1 3766 A for instance). Their 
non-detection points towards a different source than thermal an- 
nealing, which in turn means that the maximum temperature we 
find is not appropriate and that grains should be colder than 
~ 1500 K. Since the derived geometry of the disk is dominated 
by the fitting of emission features and given that there is only 
one strong emission feature around 9 //m in the entire spectrum 
of HD 15407 A, this can explain why the temperature determina- 
tion may not be accurate enough in our best model. Additionally, 
we discuss previously the possibility of a cold, outer dust belt 
further away from the central star. If confirmed, the geometry of 
the disk would have to be revised, and hence at the same time 
the maximum temperature of the grains. This temperature issue 
can be easily constrained with high-angular resolution observa- 
tions. From the absence of crystalline olivine grains, ruling out 
thermal annealing as a source of crystallization process leaves 
us with little options to explain the peculiar dust composition of 
this dust belt. Since such crystalline grains are expected to be 
common in the outer regions of primordial disks where comets 
are formed, and since silica is not usually detected in cometary 
bodies, the scenario of an outer belt feeding the inner regions is 
rather unlikely. We may therefore be witnessing the aftermath 
of collisions between silica-rich parent bodies, or the production 
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of coesite via collisions that provide temperatures and pressures 
high enough to form this particular polymorph of Si02. In the 
first case, collisions should be either non-disruptive, or the parent 
bodies should be non-porous, otherwise shock-induced heating 
may produce other crystalline grains, as discussed previously. In 
that scenario, the collision may be the consequence of dynamical 
stirring, triggered by the secondary of the binary system. 

7.5.5. Other sources 

For the three remaining sources, we can hardly discuss the origin 
of the crystalline grains or the petrology of olivine grains. For 
HD 98800 B, we conclude in Sect.[52]that the crystallinity frac- 
tion could only be considered as tentative because of the lack of 
strong emission features associated with crystalline grains. For 
HD 169666, the overall low SNR of the data renders the analy- 
sis difficult, and we prefer not to over-interpret data. Concerning 
BD-H20 307, since the fit is not satisfying around lOjum, the con- 
fidence level of the crystallinity fraction and iron content is not 
sufficient for any interpretation of their origin. 

Overall, the conclusion of this Section|7] is that many pro- 
cesses are active in these transient debris disks, from thermal an- 
nealing at high or intermediate temperatures, to (non-) disruptive 
collisions of (non-) porous parent bodies with various composi- 
tions. Spectroscopic observations contain a finite amount of in- 
formations, and we have a limited knowledge about dust emis- 
sivity properties (from near- to mid-IR). We are able to suggest 
possible scenarios to interpret our results within these unfortu- 
nate limitations but further investigations (observations, instru- 
mentation and laboratory experiments) are mandatory to con- 
firm, or infirm our findings. To summarize, it is not trivial to 
assume that the inferred dust composition originates from par- 
ent bodies inherited from the primoridal disk because of several 
other crystallization processes at stake. 

8. Comparison with Solar System objects 

In this Section, we discuss how our results compare with studies 
of bodies that were formed during the youth of our Sun. For 
instance, very few chondrites (~ 2%. [Norton & Chitwoodll2008l) 
that fall on Earth seem to contain enstatite, a result consistent 
with our modeling results. The Fe content is also an interesting 
diagnostic to focus on, as it is quite a novel finding for dust in 
exo-zodiacal belts. 

8.1. Cometary bodies 

As detailed in the reviews bv iHanner & Zolenskvl (1201 Ol) and 
IWoodenI (l2008h comets are formed in the outermost regions 
(Oort Cloud or the Kuiper Belt region) and trace the youth of 
the Solar System. Informations about their dust composition is 
a prime objective in order to better understand the formation 
of our Solar System. Using spectr oscopic observations at dif- 
ferent epochs. IWooden et al.l (1 19991) studied the spectrum of the 
comet C/1995 01 (Hale-Bopp). According to their modeling they 
found both Mg-rich olivine and crystalline pyroxene grains to be 
present in the coma, but no clear evidences for Fe-rich olivine 
grains. The Mg-rich pyroxene grains were required to reproduce 
the 9.3 jum feature, which is not detected in any of our sources. 
The modeling results for the comet Hale-Bopp a re very similar 
to the ones for the Oort Cloud comet C/2001 Q4 (IWooden et al.l 
I2OOI . For both comets, the mineralogical findings echo the re- 
sults from the Stardust mission, on comet P81/Wild2. Based 
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on samples collected in situ, IZolenskv et al.l (l2008l) found that 
olivines in the comet mostly consist of Mg-rich olivine grains, 
with a Mg fraction (=Mg/[Mg + Fe]) of about 99. However, 
they also found the composition range of olivine to be extremely 
broad, ranging from Mg fractions between 4 and 100. From the 
peak position of the Mg distribution, the author conclude that 
composition changes due to capture heating must have been in- 
significant, otherwise the Mg distribution should instead peak 
around 40-60. Interestingly, from the peak position of the Mg 
distribution, the authors suggest that the petrology of the olivine 
grains can be explained by aqueous alteration; Fe-rich olivine 
grains m ay not survive as easily as Mg-rich grains in presence 
of water (iWogelius & Walthedll992l) . The olivine grains may ei- 
ther trace the eventual presence of liquid water on P8 1 /Wild 2, or 
trace different regions of the initial disk, where Mg-rich olivine 
grains have preferentially survived aqueous alteration compared 
to Fe-rich olivine grains. The fact that crystalline grains detected 
in cometary bodies are Mg-rich rather than Fe-rich is consistent 
with observations of gas-rich Class II disks. Comets are expected 
to trace the mineralogy of the outermost regions of the primor- 
dial disk, therefore if Fe-rich olivine grains, as we detect them 
in debris disks, were already formed during the Class II phase, 
we should in principle observe significant fractions of them in 
comets as well. For the couple of objects for which we find more 
Fe-rich than Mg-rich olivine grains, this points towards a new 
generation of dust in the transient phase of the disks, which is 
subjected to alternative crystallization processes, more inclined 
to form Fe-rich crystalline grains. 

8.2. Interplanetary Dust Particles 

The exact origin of Interplanetary dust Particles (IDPs) is diffi- 
cult to track, they could have been trapped in asteroids or short- 
perio d comets (originating from the Kuiper Belt, see iBradlevI 
|20l3 for a review). It is in principle possible to distinguish be- 
tween an asteroid or comet origin, bas ed on He rele ase dur- 
ing their entry in the atmosphere (Brownl ee et al.l[T995l) . but the 
place of formation of these IDPs is thought to be presolar, based 
on isotopic considerations. An interesting point with respect to 
our study is that olivine grains contained in several IDPs ar e 
found to be Fe-rich (e.g.. lBradlevlll9. 94', 'Bru netto et all 1201 ll) . 
According to lZolenskv & Barrett! (1 1994i) the Mg distribution for 
hydrous IDPs ranges between 76 and 100, while it ranges be- 
tween 52 and 100 for anhydrous IDPs, these differenc es being 
the po ssible consequence of aqueous alteration. Ngu ven et al.l 
(I2007h discuss the origin of Fe in IDPs and suggest that the Fe 
enrichment may be the consequence of a secondary process tak- 
ing place in the solar nebu la or on the parent b ody, which echo 
our findings. Nevertheless. iNguven et al.l(l2007l) also discuss the 
possibility of condensation under non-equilibrium conditions in 
stellar outflows to explain the detection of Fe-rich crystalline 
grains. Overall, the existence of Fe-rich olivine grains in various 
astronomical environments is not to be questioned, the novelty of 
our study is to detect such grains in several exo-zodiacal belts. 
Such grains are not expected to be formed in gas-rich Class II 
disks, their detection in debris disks points toward alternative 
crystallization process during the transient phase of the disks. 



8.3. S-type asteroids 

Investigating these alternative crystallization mechanisms can be 
supported by our current knowle dge of other bodies fro m the 
Solar System, such as asteroids. iNakamura et al.l (1201 ll) study 



the mineralogy of samples collected by the Hayabusa spacecraft, 
directly on the near-Earth Itokawa S-type asteroid. The authors 
find the most abundant crystalline silicates to be of olivine com- 
position, with an averaged Fe fraction of about 30 % (Fa29.o±o.7)- 
Based on estimations of equilibration temperatures, they con- 
clude that these particles have experienced a peak temperature 
of about ~ HOOK before cooling down slowly to ~900K. To 
explain the slow cooling, the authors suggest that such particles 
shoud have been formed deep inside a large (>20km) asteroid, 
where intense thermal metamorphism can eflic iently take place 
suppo rted by the decay of Al^^ for instance). INakamura et al.l 
1201 ll) conclude that results from the Hayabusa mission pro- 
vide a link between the composition of S-type asteroids and or- 
dinary chondrites. Together with their findings, we are able to 
connect the mineralogy and petrology of bodies from our Solar 
System with dust grains observed in exo-zodiacal belts around 
Solar analogs. 

9. Conclusion 

In this study, we present a new and powerful radiative transfer 
code, dedicated to debris disks (Debra). Combining this code 
with recent laboratory experiments on dust optical properties 
we model the Spitzer/iRS spectra of 7 debris disks around solar- 
analogs, in order to derive the mineralogy of the dust. The entire 
SEDs, from optical to millimeter wavelengths, are successfully 
modeled within the optically thin regime. Based on the modeling 
of these rare objects, we obtain the following results: 

1. To contribute significantly in the mid-IR, grains must be 
"dirty", especially grains with low absorption efficiencies 
in the near-IR, such as Fe-free amorphous grains, enstatite 
and crystalline olivine grains. If Qabs values are too low 
at the wavelengths where the stellar radiation peaks, grains 
will be almost transparent and consequently too cold for 
their emission to be significant. Even though the crys- 
talline olivine grains we use in this study contain some 
Fe (M g7xFe?fi_x^Si04, with x =0.925 or 0.8). IZeidler et al.l 
( 1201 lb have demonstrated their absorption efficiencies to be 
very low in the near-IR. Adding impurities, such as carbon, 
helps circumventing this issue. One should note that we do 
not assume the dust species to be in thermal contact. Each 
grain, for all sizes and all compositions has its temperature 
determined as a function of its distance r and absorption ef- 
ficiencies. 

2. For several sources, especially HD 11 3766 A and HD 69830, 
we find that crystalline olivine grains contain Fe inclusions, 
up to 20% compared to Mg. Fe-rich crystalline grains are 
usually not observed in Class II protoplanetary disks, which 
can be expl ained by several studies of crystallization pro- 
cesses (e.g.. lGailll2004llNuth & Johnsonll2006l iMurata et"al1 
2009). Therefore our findings of more Fe-rich compared 
to Mg-rich olivine grains points towards a new generation 
of crystalline dust grains, for which the crystallization pro- 
cesses are slightly different in debris disks than in Class II 
disks (e.g., inside differentiated bodies), leading to a higher 
fraction of Fe. 

3. In the formalism of IWvatt et al.l (l2007h . using the ratio be- 
tween observed and maximal fractional luminosities, we find 
that most, if not all, of the debris disks we study are in a tran- 
sient phase, that cannot be explained by a steady-state evo- 
lution from the Class II phase, meaning that yum-sized dust 
grains were produced relatively recently compared to the 
stellar ages. Computing the radiation pressure efficiency fS^p, 
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we find that /zm-sized grains are short-lived around at least 
three objects (HD 15407 A, HD 113766 A and HD 169666). 
Investigating their variability in time via archive mining, we 
conclude there must be an on-going process that replenishes 
efficiently the population of smallest grains on short dynami- 
cal timescales. At least 5 out of 7 sources are either in binary 
systems or have planets orbiting the central star, suggesting 
the transient phase might have been triggered by dynamical 
stirring from these massive bodies. 

4. We find that relative abundances of enstatite grains are over- 
all small for all the sources we study. This is the conse- 
quence of non-detections of recognizable emission features 
associated with pyroxene crystalline grains (especially at 
9.3 fim). If the observed crystallinity is the consequence of 
present day crystallization processes, this may be explained 
by higher activation energies for annealing compared to crys- 
talline olivine grains. If the observed crystallinity fraction is 
inherited from parent bodies that collided (with no shock- 
induced heating), this echoes the scarcity of enstatite chon- 
drites in our own Solar System. 

5. The quality of the fits we present in this study is intimately 
connected to the gabs values we use in the modeling, es- 
pecially for olivine grains (which are responsible for the 
strongest emission features in the spectra we model). This 
underlines the extreme importance of laboratory experiments 
to improve our understanding of cosmic dust properties. The 
aerosol technique is a promising technique that provide di- 
rect measurements of optical properties , with as less environ - 
mental contamination as possible (e.g.. lTamanai et an2009l) . 

6. According to far-IR observations, only one source in our 
sample (HD 113766 A) harbors a cold, outer belt. Based ei- 
ther on disk properties from our modeling results (rout), or 
on mineralogical considerations (fit to emission features), 
we discuss the possibility of outer belts around two other 
sources (HD 15407 A and BD+20 307). 

As a conclusion, we would like to stress the fact that SED 
modeling is a challenging problem, for proto-planetary and de- 
bris disks. Since debris disks are optically thin at all wave- 
lengths, the problem becomes even more difficult as several as- 
sumptions no longer hold, such as thermal contact between dust 
grains (which renders the temperature determination critical) or 
the continuum emission as commonly used when modeling spec- 
tra of Class II disks. Spectroscopic observations do not contain 
any spatial informations, which introduces uncertainties in the 
modeling as demonstrated by the bayesian inference performed 
on the disk parameters. Therefore more refined modeling would 
require high angular resolution observ ations (e.g., near-IR in- 
terferometry or scattered light imaging, iLebreton et al.lEoTll) to 
lower the number of degenerate parameters (rin, rout and possi- 
bly a). Further work is also required on the dust properties to 
improve the modeling results, not only with laboratory experi- 
ments (e.g., aerosol measurements for other dust species), but 
also to investigate additional effects such as grains porosity, or 
inhomogeneous ag gregates of diff'erent species (as d escribed in 
iMin et al.l20"08l and ! Voshchinnikov & Hennind200^ . This latter 
point may eventually provide a suitable solution for disks where 
jum-sized grains are subjected to efficient radiation pressure. If 
small grains are embedded in larger grains, they may survive 
longer in the dust belt. Unfortunately, this increases the com- 
plexity of the modeling, and again underlines the need for addi- 
tional constraints on the disk geometry. With observations and 
laboratory data available at present days, we are able to success- 
fully model the Irs spectra of several warm debris disks, with 



a limited number of dust species, and provide new insights into 
the origin of this transient dust. As exemplified by Figure|9] we 
are confident in our mineralogical results, especially concerning 
to the relative abundances of Fe-rich and Mg-rich olivine grains. 

Acknowledgements. The authors thank the anonymous referee for the construc- 
tive and positive comments, as well as the A&A Editor, Malcolm Walmsley, 
for his careful reading and for providing additional suggestions to improve the 
paper J. O. acknowledges the Alexander von Humboldt foundation for critical fi- 
nancial support. A. M. and P. A. appreciate financial support from the Hungarian 
Research Fund OTKA 101393. 



References 

Abraham, P., Juhasz, A., Dullemond, C. P, et al. 2009, Nature, 459, 224 
Aumann, H. H., Beichman, C. A., Gillett, F. C, et al. 1984, ApJ, 278, L23 
Begemann, B., Dorschner, J., Henning, T., Mutschke, H., & Thainin, E. 1994, 
ApJ, 423, L71 

Beichman, C. A., Bryden, G., Gautier, T. N., et al. 2005, ApJ, 626, 1061 
Beichman, C. A., Lisse, C. M., Tanner, A. M., et al. 2011, ArXiv e-prints 
Bouwman, J., Henning, T, Hillenbrand. L. A., et al. 2008, ApJ, 683, 479 
Bouwman, J., Meeus, G., de Koter, A., et al. 2001, A&A, 375, 950 
Bradley, J. 2010, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 815, 

Lecture Notes in Physics, Berlin Springer Verlag, ed. T. Henning, 259-276 
Bradley, J. P 1994, Science, 265, 925 

Britt, D. T., Yeomans, D., Housen, K., & Consolmagno, G. 2002, Asteroids III, 
485 

Brownlee, D. E., Joswiak, D. J., Schlutter, D. J., et al. 1995, in Lunar 
and Planetary Institute Science Conference Abstracts, Vol. 26, Lunar and 
Planetary Institute Science Conference Abstracts, 183 
Brunetto, R., Borg, J.. Dartois, E., et al. 201 1, Icarus, 212, 896 
Bryden, G., Beichman, C. A., Trilling, D. E., et al. 2006, ApJ, 636, 1098 
Bums, J. A., Lamy, P L., & Soter, S. 1979, Icarus, 40, 1 
Carpenter, J. M., Bouwman, J., Mamajek, E. E., et al. 2009, ApJS, 181, 197 
Charbonneau, R 1995, ApJS, 101, 309 

Chen, C. H., Mamajek, E. E., Bitner, M. A., et al. 2011, ApJ, 738, 122 
Chen, C. H., Patten, B. M., Werner, M. W., et al. 2005, ApJ, 634, 1372 
Chen, C. H., Sargent, B. A., Bohac, C, et al. 2006, ApJS, 166, 351 
Davison, T. M., Collins, G. S., & Ciesla, R J. 2010, Icarus, 208, 468 
Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531 

Dorschner, J., Begemann, B., Henning, T., Jager, C, & Mutschke, H. 1995, 

A&A, 300, 503 
Egan, W. G. & Hilgeman, T. 1977, Icarus, 30, 413 

Fabian, D., Jager, C, Henning, T., Dorschner, J., & Mutschke, H. 2000, A&A, 
364, 282 

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2012, ArXiv e- 
prints 

Fujiwara, H., Onaka, T, Ishihara, D., et al. 2010, ApJ, 714, L152 
Furlan, E., Sargent, B., Calvet, N., et al. 2007, ApJ, 664, 1176 
Gail, H.-P 2004, A&A, 413, 571 

Gail, H.-P. 2010, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 815, 

Lecture Notes in Physics, Berlin Springer Verlag, ed. T. Henning, 61-141 
Gorlova, N., Balog, Z., Rieke, G. H., et al. 2007, ApJ, 670, 516 
Hanner, M. S. & Zolensky, M. E. 2010, in Lecture Notes in Physics, Berlin 

Springer Verlag, Vol. 815, Lecture Notes in Physics, Berlin Springer Verlag, 

ed. T. Henning, 203-232 
Henning, T. 2010, ARA&A, 48, 21 
Henning, T. & Mutschke, H. 1997, A&A, 327, 743 
Henning, T. & Mutschke, H. 2010, Journal of Nanophotonics, 4, 041580 
Jager, C, Dorschner, J., Mutschke, H., Posch, T., & Henning, T. 2003, A&A, 

408, 193 

Jager, C, Molster, F. J., Dorschner, J., et al. 1998a, A&A, 339, 904 

Jager, C, Mutschke, H., & Henning, T. 1998b, A&A, 332, 291 

Juhasz, A., Bouwman, J., Henning, T, et al. 2010, ApJ, 721, 431 

Juhasz, A., Dullemond, C, van Boekel, R., et al. 2011, ArXiv e-prints 

Keil, K., Stoeffler, D., Love, S. G., & Scott, E. R. D. 1997, Meteoritics and 

Planetary Science, 32, 349 
Kemper, F, Vriend, W. J., & Tielens, A. G. G. M. 2005, ApJ, 633, 534 
Kitamura, R., Pilon, L., & Jonasz, M. 2007, Appl. Opt., 46, 8118 
Koemer, D. W., Jensen, E. L. N., Cruz, K. L., Guild, T. B., & Gultekin, K. 2000, 

ApJ, 533, L37 

Koike, C, Chihara, H., Tsuchiyama, A., et al. 2003, A&A, 399, 1 101 
Koike, C, Imai, Y, Chihara, H., et al. 2010, ApJ, 709, 983 
Krivov, A. V. 2010, Research in Astronomy and Astrophysics, 10, 383 
Lebreton, J., Augereau, J.-C, Thi, W.-F., et al. 201 1, ArXiv e-piints 
Liseau, R., Eiroa, C, Fedele, D., et al. 2010, A&A, 518, L132 



24 



Olofsson et al.: Transient dust in warm debris disks 



HD113766A 



HD69830 





2.5 




z . u 




1.5 












1.0 




0.5 




o.oL 








40- 




20 







b " 


-20 




-40 




10 



10 



15 



20 



25 



15 



20 



25 



BD+20 307 



20 25 

A [^m] 



30 



30 



35 



35 




0.15 



^ 0.10 



0.05 



O.OO 




„ 40- 
^ 20 

■q -20 

-40- 



O.OL 



20 
10 

0-\ 
-10 
-20 




10 



10 



10 



15 



20 25 

A [^m] 

HD15407A 



30 



15 



20 



25 



30 



15 



20 25 

A Ifim] 

HD98800B 



30 



35 




35 



35 



40 



40 



40 




60 
T 40 

5S 20 

r 

1-20 
° -40 

-6a 



10 



IDS 




15 



20 

A [^m] 



25 



30 



35 



25 



Fig. 14. Best fits (in red) to the Irs spectra (dashed lines). Uncertainties are 3-cr. For each source bottom panels show the residuals. 
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Fig. 15. Complete SEDs for all the sources, with Irs spectra and photometric observations in red. Photospheres are in dashed black 
and modeled SEDs are in solid black. 
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